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ABSTRACT 


The  inherent  symmetry  of  a  spherical  reflector  microwave  antenna 
suggests  that  if  the  fields  incident  upon  the  reflector  from  a  feed 
required  to  produce  a  specified  reflected  field  can  be  determined, 
then  scanning  of  this  reflected  field  may  be  accomplished  by  rotation 
of  the  feed  assembly  rather  than  by  rotation  of  the  entire  antenna. 
A  technique  has  been  determined  to  solve  for  the  fields  in  vicinity 
of  a  hypothetical  transmitting  feed  that  will,  following  reflection 
from  the.  surface  of  the  spherical  reflector,  yield  a  specified  field 
pattern  at  the  reflector  aperture.  From  given  fields  over  a  sphere, 
a  portion  of  which  coincides  with  the  surface  of  the  spherical  reflec- 
tor, a  spherical  harmonic  series  expansion  is  used  to  calculate  the 
fields  on  a  smaller  concentric  sphere  encompassing  the  feed  region 
in  solution  of  the  boundary  value  problem. 

This  technique  is  extended  to  a  case  where  it  is  shown  that 
geometrical  optics  approximations  are  valid  over  a  specific  region 
of  the  reflector  system,  and  that  by  ray  tracing,  the  fields  on  an 
arbitrary  reference  sphere,  whose  radius  is  confined  only  to  the 
axis  of  symmetry  of  the  reflector,  can  be  determined.  From  this 
field  distribution,  the  spherical  harmonic  series  expansion  can  be 
used  to  determine  the  fields  on  the  surface  of  a  smaller  concentric 
sphere  enclosing  the  hypothetical  feed.  It  is  shown  that'  this  tech- 
nique has  flexibility  wherein  a  number  of  origins  could  be  examined 
for  various  feed  spheres  with  fewer  terms  of  the  series  expansion 
than  with  the  fixed-origin  case. 
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CHAPTER  I 
GENERAL   BACKGROUND 

1.1    The  fact  that  a  reflecting  paraboloid  of  revolution,  when  il- 
luminated by  a  source  of  spherical  waves  placed  at  the  geometrical 
focus  of  the  paraboloid,  will  produce  a  uniform  plane  wave  over  the 
aperture  of  the  reflector  is  veil  known.  The  generation  of  a  pencil- 
shaped  beam  antenna  pattern  has  reached  a  state  of  refinement  such 
that  extremely  narrow  bcamwidth,  lev;  sidelobe  levels,  reasonably 
broad-band  operation,  high  aperture  efficiency  and  gain  are  readily 
attainable  for  many  applications. 

The  axis  of  this  pencil  beam  coincides  with  the  axis  of  revolution 
of  the  reflector,  and  attempts  to  direct  the  antenna  pattern  by  off- 
axis  movement  of  the  point-source  feed  have  failed  for  all  but  the 
smallest  off-axis  angles.  Therefore,  to  direct  the  pattern  requires 
that  the  entire  reflector  and  feed  assembly  be  rotated  or  moved  to  aim 
the  antenna  pattern  in  the  direction  of  interest. 

Widespread  use  of  large  aperture  paraboloidal  reflectors  with  wide- 
angle  scanning  requirements  has  stimulated  considerable  interest  in 
development  of  improved  feed  systems  to  alleviate  problems  of  ground 
noise  radiation,  spurious  polarization  effects  and  low  aperture  ef- 
ficiencies encountered  in  radio-telescope  applications.  Other  such 
problems  that  are  encountered  include  difficulties  in  maintaining 
mechanical  and  structural  tolerances  of  the  reflecor  surface  at  partic- 
ular frequencies  of  interest. 


The  inherent  symmetry  of  the  spherical  reflector  suggests  that  it 
could  provide  some  relief  from  the  problems  previously  described.  Assum- 
ing an  energy  source  could  be  obtained  that  would,  upon  reflection  from 
the  spherical  surface f  provide  an  antenna  pattern  whose  beam  coincided 
with  the  radius  of  the  reflector  upon  which  illumination  occurred,  it 
would  be  no  longer  necessary  to  overcome  the  inertia  of  the  entire 
system  to  meet  a  wide-angle  scanning  requirement,  but  merely  to  direct 
the  feed  assembly  as  desired.  In  such  a  case,  it  is  feasible  that  large 
apertures  could  be  illuminated  to  provide  desired  gain  .  Further,  the 
symmetry  of  a  spherical  section  enhances  both  ease  of  fabrication  and 
verification  of  system  tolerances  in  many  applications.  Other  advantages 
that  may  be  realized  from  such  a  system  follow  directly  from  speculation 
of  possible  application,  including  for  example  an  antenna  system  whose 
reflector  could  be  imbedded  in  the  earth,  with  few  moving  parts,  pre- 
senting a  low  silhouette  to  wind  forces. 

Unfortunatly,  analysis  of  the  spherical  reflector  reveals  its 

inherent  spherical  aberration.  For  a  plane  wave  incident  upon  a 

spherical  reflector,  rays  close  to  the  axis  of  symmetry  intersect  at 

a  point  defined  as  the  paraxial  focus,  exactly  half-way  between  the 

center  of  curvature  of  the  spherical  section  and  the  vertex.  Other  rays 

reflected  from  the  surface  intersect  the  axis  at  points  closer  to  the 

vertex.  This  phenomenon  may  be  observed  in  figure  1,  where  rays  are 

constrained  by  Snell's  law  to  reflect  with  an  angle  equal  to  that  of 

incidence,  as  measured  from  the  surface  normal  at  the  point  of  interest. 

A  more  complete  explanation  of  this  reflection  phenomenon  will  be  given, 
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This  suggests  that  rather  than  a  distinct  point  focus  there  exists  a 
fuzzy  focal  region.  Illumination  of  the  reflector  with  a  point  source  in 
the  focal  region  reveals  a  diffraction  pattern  with  wider  beamwidth, 
higher  sidelobes  and  much  less  gain  than  a  paraboloidal  reflector  of 
identical  aperture.  It  is  this  aberration  that  has  presented  an  obstacle 
to  development  of  the  spherical  rather  than  paraboloidal  reflector  for 
use  as  a  microwave  antenna?  Since  introduction  of  the  problem  to  the 
microwave  cortmunity  in  World  War  II,  however,  there  has  remained  an 
active  interest  in  the  problem.  Most  recently,  with  construction  in 
1960  of  the  Arccibo  (Puerto  Rico)  Ionospheric  Observatory's  1000-foot 
diameter  spherical  reflector  ,  considerable  interest  has  been  generated 
in  the  solution  to  problems  of  spherical  aberration  at  radio  frequencies 
by  means  of  various  feed  systems. 
1,2   It  is  appropriate  to  review  past  efforts  relevant  to  the  problem.  Early 
attempts  to  explain  optical  phenomena  followed  from  geometrical  theory. 
Euclid  presented  an  analysis  cf  the  reflection  laws  around  300  B.C.  while 
observing  that  light  energy  could  be  concentrated  by  means  cf 


2 
reflectors.  Early  astronomers  used  spherical  mirrors  as  primary 

collectors  of  light  incident  upon  their  telescopes.  Presumably  the 
ease  of  manufacture  of  these  spherical  mirrors  led  to  their  V;idespread 
use,  even  though  the  geometrical  analysis  of  the  reflecting  properties 
of  paraboloidal  surfaces  had  been  conducted.  The  difficulty  of  manufacture 
of  the  latter  apparently  precluded  exploiting  their  point-focus  charac- 
teristic. 

The  concept  of  a  spherical  reflector  for  microwave  use  evolved  from 
the  World  War  II  advances  of  radar  technology.  In  1941,  Ashmead  and 
Pippard  suggested  two  approaches  to  the  aberration  encountered: 

1.  For  reflectors  of  large  radius  and  restricted 
aperture,  the  surface  of  the  spherical  reflec- 
tor departs  only  slightly  from  that  of  a  para- 
boloidal reflector.  It  follows  that  the  diffrac- 
tion pattern  is  only  slightly  different  from  that 
of  a  paraboloidal  reflector.  If  the  slight  depar- 
ture from  plane  wave  illumination  at  the  aperture 
and  the  value  of  gain  accompanying  the  restric- 
tive aperture  can  be  tolerated,  this  approach 
permits  a  simple  design  and  is  suitable  for  wide 
angles  of  scan. 


2,  For  the  general  case  of  spherical  reflectors „ 
should  the  primary  energy  source  be  a  special 
radiator  that  provides  not  spherical  waves  but 
waves  whose  departure  from  spherical  waves  is 
just  that  necessary  to  compensate  for  the  differ- 
ence between  the  sphere  and  paraboloid,  then 
reflected  wave  fronts  will  be  plane  waves,  uniform- 
ly illuminating  the  aperture. 

Early  work  indicated  that  it  was  possible  to  create  non-spherical 
waveo  of  approximately  the  right  shape  by  proper  combination  of  wave- 
guide horns  and  dielectric  lenses.  Unfortunatly,  difficulties  in 
Correcting  for  the  larger  phase  error  did  not,  at  the  time,  justify 
pursuance  of  this  approach.  Rather,  to  meet  the  expedience  of  the 
wartime  effort,  the  first  approach  of  toleration  of  error  rather  than 
compensation  was  used  exclusively.  Subsequent  studies  exploiting  this 
technique  reached  a  natural  conclusion  with  the  work  of  Li  ,  who 
extended  and  verified  the  theory  of  restrictive  aperture  large  radius 
spherical  reflectors.  His  experimental  work  indicated  that  with  a  ten- 
foot  aperture  hemisphere  absolute  gain  comparable  to  that  of  a  forty- 
inch  paraboloid  was  obtained  over  a  scan  of  140  with  rotation  of 
the  feed  at  a  constant  radius  from  the  reflector. 

The  second,  more  interesting  approach  to  the  solution  allows  the 
use  of  a  larger  aperture  at  the  expense  of  increased  complexity  in 

cost,  design  and  some  reduction  in  the  available  angle  of  scan. 
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In  the  context  of  the  latter  approach  corrective  devices  have  fallen 
into  the  following  categories: 

1.  lens  structures 

2.  auxilary  reflectors 

3.  line-source  feeds 

4.  transverse  aperture  feeds 

Lens  structures,  while  providing  satisfactory  correction  at  optical 
frequencies,  provide  marginal  correction  in  microwave  applications.  Sizes 
comparable  to  the  reflector  and  primary  source,  requirements  for  large 
pieces  of  carefully  controlled  precision-cut  dielectric  and  inherent 
dielectric  losses,  to  say  nothing  of  the  costs  and  weights  involved, 
relegate  this  category  of  corrective  devices  at  microwave  frequencies  to 
a  position  of  academic  interest, 

Auxilary  reflectors  are  of  minor  interest  when  the  undesirable  charac- 
teristics of  shadowing  or  aperture  blockage  are  imminent.  Even  in  cases 
where  this  is  of  small  concern,  the  excessive  weight  ef  the  secondary 
reflector  must  be  considered,  as  must  the  mechanical  arrangements  for 
support  of  the  assembly,  A  gregorian  reflector  has  been  described  by 
Holt  and  Bouche  wherein  a  reflecting  surface  intercepting  a  great  deal 
of  the  energy  incident  from  the  spherical  reflector  is  installed.  The 
surface  of  the  sub-reflector  is  such  that  it  provides  a  constant  path 
length  from  seme  reference  plane  near  the  aperture  to  the  primary  source. 


Experimental  results  indicate  that  the  scanning  capability  was  realized 
to  30  by  rotation  of  the  sub-reflector  on  a  constant  radius  from  the 
spherical  surface,  Sidelobe  levels  were  very  high  in  comparison  to 
those  available  with  paraboloidal  reflectors „  While  aperture  blockage 
was  estimated  to  be  1,5%,  no  figures  on  gain  were  given. 

The  earliest  attempts  at  compensation  involved  the  fact  that  rays 
incident  upon  the  spherical  reflector  were  reflected  to  cross  the  axis 
of  symmetry.  It  was  postulated  by  Spencer,  Sletten  and  Walsh  that  the 
incoming  energy  could  be  recombined  with  an  appropriate  phase  shift 
to  provide  a  maximum  of  received  signal.  Reciprocity  was  invoked  to 
suggest  that  a  properly  phased  line  source  could  produce  a  wavefront 
that  would  provide  uniform  illumination  over  the  aperture.  Techniques 
of  geometrical  optics  were  used  for  path  length  calculations.  The 
sources  considered  were  open-ended  waveguides,  polyrod  travelling  wave 
structures,  discrete  dipole  sources  alog  the  axis  and  slotted  waveguide 
feeds.  Indeed t  the  first  feed  system  constructed  at  the  Arecibo 
reflector  was  based  upon  a  modification  of  this  design.  Aberrationless 
scan  of  over  60  was  predicted  for  this  initial  installation-as  well 
as  reasonable  aperture  efficiency.  Actual  results  were  veil  under  those 
anticipated.  The  high-power  dual  circularly  polarized  feed  illuminated 
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the  full  aperture  of  1000  feet  with  an  efficiency  of  21%  at  430  MHz. 

9 
Love  described  a  composite  radiating  system  wherein  channel  guide 

feed  and  a  polyrod  array  provided  a  travelling  wave  structure  to  illumin- 
ate from  both  the  line  focus  and  the  paraxial  focus.  He  predicted  a  scan 
of  110  for  a  hemispherical  reflector,  and  efficiency  on  the  order  of 

12. 


55%«  Schell  investigated  the  fields  along  the  axis  by  application  of 
diffraction  theory.  He  demonstrated  the  deviation  from  the  geometrical 
optics  solution  in  the  non-zero  wavelength  and  agreement  in  the  zero- 
wavelength  case.  Schell1 s  solution  consisted  of  the  use  of  a  sufficient 
number  of  discrete  sources  located  along  the  axis,  illuminating  partic- 
ular regions  of  the  reflector.  The  analysis  was  based  upon  the  fields 
present  along  the  focal  axis  when  the  reflector  was  illuminated  by  an 
incoming  plane  wave.   In  an  attempt  to  solve  the  problem  based  upon 
Schell' s  solution,  Cohen  and  Perona   constructed  a  feed  system  at  the 
Arecibo  reflector.  Their  compound  antenna  feed  consisted  of  a  four- 
element  linear  array  with  a  reflecting  screen  and  a  sixteen-element 
linear  array  located  on  the  axis.  Experimentation  determined  that  gain 

was  lower  and  sidelobes  higher  than  had  been  anticipated.  Aperture 
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efficiencies  were  on  the  order  of  30-38%   .  McCormick   followed  this 

study  with  an  analysis  that  considered  the  transmitting  case  and,  in 
particular,  the.  aspect  of  polarization.  His  analysis  concluded  that  a 
gain  deterioration  was  inevitable  with  a  circular  waveguide  with  cir- 
cumferential slots  that  did  not  consider  the  azimuthal  fields  as  well 
as  the  longitudinal  fields.  The  theory  of  McCormick  was  used  by 
Love  and  Gustinic   to  illustrate  the  feasibility  of  illuminating  from 
a  leaky  cylindrical  waveguide  that  considered  these  fields.  Host 

recent  literature  indicates  that  a  line  source  feed  with  high  performance 

15 

has  finally  been  achieved  at  Arecibo.  LaLonde  and  Harris  "  note  that 

their  feed  has  been  installed  at  the  observatory  and  has  realized  an 
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overall  efficiency  of  70%  at  318  MHz,  Designing  for  a  flat  waveguide  line 
feed  vith  transverse  slots  was  necessary  to  avoid  problems  of  previous 
feeds  which,  it  was  felt,  did  not  provide  sufficient  control  over 
amplitude  and  phase  variation  along  the  line  focus  nor  sufficient 
isolation  from  the  effects  of  mutual  coupling  between  the  radiating 
elements  of  the  feed.  The  lower  frequency  of  operation  was  specified  to 
avoid  losses  due   to  surface  imperfections  of  the  reflector.  The  design 
of  the  phase  and  amplitude  required  was  based  upon  a  geometrical  optics 
analysis  of  the  rays  intersecting  the  axis.  This  design  is  significant 
in  that  it  represents  the  first  installation  of  a  corrective  feed  in 
a  spherical  reflector  that  has  yielded  aperture  efficiencies  approaching 
these  of  paraboloidal  reflectors.  Aside  from  strengthening  arguments  for 
fixed-ref lector  scenning  antennas,  providing  some  40  of  scan  with  high 
efficiency  at  Arecibo,  the  success  of  this  feed  indicates  the  potential 
of  the  spherical  reflector  antenna  for  fully  steerable  arrays,  where 
desired. 

It  has  been  noted  by  Ricardi   that  line  source  radiator  correcting 
feed  systems  must  be  capable  of  producing  a  wave  whose  polarisation  and 
shape  arc  not  easily  obtained  with  elementary  radiators. This  basic  fact 
has  caused  present  day  line  source  correcting  feeds  to  fall  short  of 
expectations,  with  exception  of  the  LaLonde- Harris  feed.  It  should  be 
noted,  however,  that  the  requirements  of  Arecibo  call  for  high  pox^er 
with  capability  for  dual  circular  polarization.  The  line  feed  of  LaLonde 
and  Harris  achieved  success  in  providing  for  linear  polarization.  While 
the  high  power  requirement  may  be  easily  met  with  the  slotted  guide, 


the  authors  note  that  polarization  flexibility  has  yet  to  be  developed. 
Presumably,  the  transverse  aperture  corrective  feed  designs  have  no 
inherent  difficulty  in  polarization  flexibility..  Early  attempts  to 
provide  aperture  corrective  feeds  mat  with  reasonable  success  in  over- 
coming the  aberration,  but  were  hampered  by  aperture  blockage  of  the 

18 
moderately  large  feed  systems.  Sletten  and  Mavroides   were  succesful 

in  lowering  sidelobes  substantially,  while  plagued  with  aperture 

blockage.  An  aperture  feed  was  suggested  for  the  spherical  reflector  by 

19 
Burrows  and  Ricardi   in  which  they  predicted  a  possible  gain  of  61  db. 

compared  to  the  maximum  gain  of  such  an  aperture  the  size  of  Arecibo  at 
62«8  db.  Their  feed  systera  was  derived  from  a  spherical  harmonic  expan- 
sion of  the  fields  necessary  to  produce  a  uniform  distribution  over  the 

aperture,  based  upon  the  fields  incident  upon  a  spherical  reflector. 

20 
Spencer  and  Hyde   conducted  a  detailed  study  of  the  focal  region  fields 

associated  with  the  spherical  reflector  illuminated  by  incoming  plane 
waves  and  concluded  that  only  a  vector  field  solution  considering  polar- 
ization would  yield  adequate  information.  They  solved  the  diffraction 
integral  by  means  of  a  stationary  phase  approximation,  noting  the 
similarities  to  the  geometrical  optics  approximations.  Their  theory 
was  shown  to  have  reasonable  agreement  with  fields  actually  found  in  the 

focal  region  by  experimentation,  and  their  work  has  been  used  as  a 

21 

basis  for  further  study  into  aperture  feeds  .  Their  work,  however,  did 

not  offer  a  solution  to  the  problem,  but  rather  provided  a  theoretical 

analysis  of  the  fields  that  could  be  expected.  Work  by  Minuet  and 

22 

MacA.  Thomas   has  been  extended  to  show  that  the  fields  in  the  image 

is 


of  symmetrical  focussing  reflectors  can  be  generated  by  certain  hybrid 

23 

modes  to  provide  theoretically  aberrationless  behavior "  .  An  aperture 

source  is  described  that  uses  a  corrugated  waveguide1,  structure  to  genei 
a  set  of  waves  that  combines  to  provide  a  theoretically  high  efficiency. 
Further  developing  his  earlier  ideas  of  spherical  harmonic  expansions, 
Ricardi  has  recently  provided  a  technique  for  synthesis  of  the  fields 
on  a  surface  enclosing  a  transmitting  feed.  This  treatment  for  the 
transmitting  feed  synthesis  provides  for  the  solution  to  an  electro- 
magnetic boundary  value  problem,  identifying  as  boundaries  the  reflector 
surface  and  the  surface  of  a  smaller  sphere  enclosing  the  feed  region. 
It  is  with  this  technique,  substantially  different  from  those  involving 
geometrical  optics  or  approximations  of  stationary  phase,  offering  a 
solution  that  appears  to  contain  inherent  accuracy,  with  which  this 
thesis  is  concerned. 
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CHAPTER  II 
PROBLEM  STATEMENT  AND  THE  THEORY  OF  GEOMETRICAL  OPTICS 

2.1  The  solution  proposed  by  Ricardi  to  the  problem  of  the  spherical 
reflector  aicrowave  antenna  relates  two  sets  of  electromagnetic  fields 
over  portions  of  surfaces  of  two  concentric  spheres.  Due  to  the  mathe- 
matical formulation  of  the  problem,  one  surface  coincides  with  that  of 
the  reflector  while  the  other  surface  is  taken  at  a  radius  of  interest 
within  that  of  the  reflector.  Origin  of  both  spheres  is  at  the  cenLer 
of  curvature  of  the  reflector  (figure  2).  The  fields  over  the  reflector 
surface  are  specified  to  be  those  fields  that  will  produce  a  desired 
field  distribution  over  the  aperture  of  the  reflector  follcwing  reflection 
from  its  surface.  The  fields  at  the  reflector  surface  arc  then  expanded 
in  a  spherical  wave  expansion.  Thus  it  nay  be  determined  that  if  the  field 
distribution  given  by  this  expansion,  evaluated  at  the  inner  surface,  can 
be  somehow  generated  by  a  source  within  this  surface  that  the  appropriate 
reflected  fields  will  appear,  as  desired,  over  the  aperture.  Interpreta- 
tion of  the  fields  at  various  radii  has  revealed  that  there  are  certain 
radii  coincident  with  nearly  constant  phase  fronts.  Synthesis  of  these 
fields  would  be  enhanced  due  to  a  potential  lessening  of  the  complexity 
of  the  radiating  structure.  Exactly  which  radius  to  be  chosen  for  the 
inner  surface  is  beyond  the  scope  of  this  work,  but  is  covered  by 
Ricardi  . 
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It  is  the  hypothesis  of  this  v/ork  that  the  technique  of  solution 
to  the  boundary-value  problem  may  be  applied  to  a  reference  sphere 
whose  otalgia  is  a  variable  distance  along  the  axis  of  symmetry  and  to 
the  surface  of  which  outgoing  rays  from  a  hypothetical  feed  located 
within  this  reference  sphere  would  be  traced  if  the  spherical  reflector 
were  absent.  The  rays  outgoing  from  the  region  of  the  hypothetical  feed 
would  be  precisely  those  rays  representing  the  fields  that  would,  upon 
reaction  with  the  surface  of  the  spherical  reflector,  reflect  a  desired 
illumination  over  the  reflector  aperture*  These  rays  would  be  based  upon 
the  rays  reflecting  from  the  spherical  surface  when  illuminated  by  an 
incident  plane  wave  over  the  aperture.  Strictly  speaking,  this  analysis 
begins  with  the  principles  of  the  receiving  case,  but  then  requires 
the  field  distribution  to  be  produced  over  the  reference  sphere  by  some 
transmitting  feed.  It  is  then  the  reaction  of  the  radiated  field  with 
the  spherical  reflector  that  theoretically  produces  a  desired  wavefront 
over  the  aperture.  The  geometry  for  this  case  is  given  by  figure  3. 

Thus,  within  the  limits  imposed  on  ray  tracing  technique  by  theory 
of  geometrical  optics,  Ricardi's  results  for  the  transmitting  feed 
case  may  be  compared  to  an  approximate  solution  that  originates  with 
the  receiving  situation.  It  is  further  possible  that  such  a  solution, 
known  to  be  strictly  valid  only  in  the  zero-wavelength  limit,  may  yield 
results  with  less  computation  time  than  the  exact  solution  within  an 
appreciably  small  deviation  in  results.  The  choice  of  origin  of  the 
system  of  concentric  spheres  suggests  that  a  more  detailed  analysis  of 
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the  fields  in  a  particular  region  nay  be  conducted  in  an  economical 
fashion,  since  the  fields  will  be  distributed  over  a  sphere  of  smaller 
radius  than  the  previous  problem's  inner  sphere  with  origin  at  that 
of  the  reflector, 

Ue  thus  proceed  to  describe  the  illumination  of  a  spherical 
reflector  by  a  uniform  plane  wave  to  illustrate  certain  character- 
istics of  the  reflector.  This  is  followed  by  development  of  geometrical 
optics  principles  essential  to  the  description  of  the  fields  over 
the  reference  sphere.  Finally,  expressions  will  be  developed  that 
yield  the  amplitude,  phase,  and  direction  of  the  electric  field 
over  the  reference  sphere. 


2.2   It  has  been  noted  previously  that  rays  reflected  from  the  spher- 
ical surface  following  illumination  by  a  uniform  plane  wave  over  the 
aperture  cross  the  axis  of  symmetry  at  varying  distances  from  the 
paraxial  focus.  This  phenomenon  may  be  analytically  formulated  to 
illustrate  the  generation  of  caustic  surfaces  by  these  reflected  rays, 
In  fact,  it  will  be  shown  that  the  caustic  surface  is  the  envelope 
of  reflected  rays.  Consider  the  situation  given  by  figure  A. 
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Constrained  only  by  the  fact  that  the  angles  of  reflection  and  incidep.ee 
are  identical,  the  equations  of  rays  may  be  given  by  the  slope  intercept 
form  of  the  straight  line  equation.  Thr.t  the  rays  are  straight  lines  aad 
the  angles  equal  will  be  covered  in  detail  in  later  sections.  Viewing 
the  system  in  the  x-z  plane  we  have: 


x  -  x 


z  -  z 


tan  20 


(2.1) 


Points  on  the  reflector  surface  are  given  by: 

x  -  R  sine 
z  =  R  cosG 


Lotting  points  on  the  envelope  be  x  ,  z  ,  and  substituting  into  2.1 

2i 


x  -   R  sinG  -  (R  cosO  -  z  )tan20  (2, 2) 

c  c 


Thus  points  on  the  caustic  surface  as  viewed  in  the  x-z  plane  become 
functions  of  parameter  0.  We  then  form  the  function  f(x,z,0)  =  0  that 
will  be  satisfied  by  points  on  the  caustic. 


f(x,z,0)  =  x  -  R  sinG  +  (R  cosG  -  z  )tan20  =  0  (2.3) 


Adjacent  points  on  the  caustic  generated  by  rays  incident  at  angle 
0  +60  will  still  satisfy  the  function  for  snail  60.  Siiice  both  functions 
f(x  ,z  ,0)  and  f(x  ,z  r0-!-50)  are  equal  to  zeros  their  difference  as  well 
is  equal  to  zero*  We  can  then  form  the  expression  valid  in  the  limit  of 
small  60: 


f(x  ,z  ,0+60)  -  f(x  ,z  ,0)     _ 
i  j     r^  r,  c*  c'  c*  c*      =0 

Iim60->O  — — — — — — 

60 


or  d  f(x  ,z  ,0) 


Differentiating  and  solving  for  z  yields: 


zc  =  R  (2  sin20  +  l)cos0  (2.5) 

2 


2.2, 


This  expression  for  z  is  then  substituted  into  2.2  to  solve  for  x  . 
r  c  c 

Manipulation  with  trigonometric  identities  yields  the  parametric  form 
of  the  equations  for  the  caustic  surface: 


x  ■  R  sinG  =  R(  3  sin0  -  sin30)  (2.6) 

C  A 


z  ■  Rcos  (  2  sin20  +  1)  -  R(  3  cos0  -  cos30)   (2.7) 
C   2  A 


Thus  it  has  been  shown  that  the  reflected  rays'  envelope  forms  the 
caustic  surface  and  quantitative  expressions  for  this  surface  have  been 
developed.  They  will  be  called  upon  in  conjunction  with  development  of 

the  expression  for  amplitude.  


2c3   Electromagnetic  field  problems  may  be  formulated  in  two  contrasting 
fashions,  following  from  a  choice  of  application  of  the  Maxwell  field 
equations.  In  one  treatment,  relationships  may  be  derived  that  determine 
what  fields  will  arise  from  a  prescribed  set  of  sources.  The  field  and 
vector  Helmholtz  equations  are  integrated  by  application  of  a  vector 
Green's  theorem,  and  express  the  fields  at  an  observation  point  as  the 
sura  of  contributions  from  the  sources  distributed  throughout  a  particular 
volume  and  from  fields  on  the  surface  of  the  volume  arising  from  sources 
outside  the  volume.   This  rigorous  treatment  indicates  that  the  field 
represents  a  flow  of  energy  outward  from  the  region  of  the  sources,  and 
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that  these  sources  and  fields  satisfy  certain  assumptions  of  continuity 
over  the  surface  of  interest,  as  required  by  application  of  the  Green's 
theorem.  In  the  second  treatment,  the  integral  relations  derived  by  the 
first  method  may  be  applied  to  yield  expressions  for  the  field  vectors 
at  a  specified  point, given  the  values  of  the  fields  over  a  surrounding 

equiphase  surface,  or  wavefront,  without  direct  reference  to  the  sources 
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generating  the  wavefronts .  These  integral  expressions  are  the  analytic 

formulation  of  the  Huygens-Fresnel  principle,  which  states  that  each  point 
on  the  equiphase  surface  can  be  regarded  as  a  secondary  source  of  elec- 
tric and  magnetic  current  and  charge  which  in  turn  gives  rise  to  a  spher- 
ical wavelet.  The  fields  at  a  point  of  interest  may  be  obtained  by  a 
superposition  of  these  wavelets,  considering  the  phase  difference  upon 
arrival  at  the  point  of  interest.  Hence,  energy  flow  is  a  wave  phenomenon- 
one  which  may  be  studied  without  reference  to  the  sources  of  energy.  As 
before,  the  rigor  of  the  solution  depends  upon  adherence  of  the  given 
field  distribution  to  several  constraints  as  imposed  by  application  of  the 
Green's  theorem. 

While  the  Huygens-Fresnel  relation  provides  expressions  for  solutions 
to  the  wave  equation  that  satisfy  Maxwell's  equations,  complexities  of 
integration  of  these  expressions  often  render  solutions  impractical,  or 
require  approximations  to  be  made  that  eventually  lessen  the  rigor  of  the 
solution  only  after  substantial  manipulation.  In  these  cases  it  becomes 
convenient  to  approach  the  propagation  of  waves  from  the  standpoint  of 
geometrical  optics,  wherein  a  rigorous  solution  is  sacrificed  for  a 


simpler  formulation  that  offcrr.  a  more  immediate  solution,  subject,  to 
several  additional  constraints  on  the  problem.  In  geometrical  optics, 
successive  positions  of  cquiphase  wavefronts  of  the  fields  and  an 
associated  system  of  rays  are  related. 


In  an  arbitrary  medium  the  wave  field  is  characterized  by  both  a  ray 
velocity  and  wave  velocity  at  every  point,  where  the  ray  velocity  is  the 
velocity  of  energy  propagation,  directed  as  the  ray  passing  through  the 
point  of  interest,  and  the  wave  velocity  is  the  rate  of  displacement  of 
the  wavefront  in  the  direction  normal  to  the  surface.  The  energy  propa- 
gates  at  the  speed  of  light  in  vacuum,  c.  The  wavefront  propagates  at  the 
velocity  v,  a  function  of  the  medium. 

We  proceed  to  describe  the  relations  among  wavefront  and  phase,  con- 
sidering vavefronts  generated  by  some  source  propagating  in  a  steady- 
state  radiation.  Let  the  wavefront  at  time  t  be  the  surface  L(xsy,z)--:L  , 

o  o* 

a  surface  of  constant  phase  V  =  (to/c)L  relative  to  some  point.  The  wave- 
front  at  some  short  time  6t  latey  is  the  surface  L(x,y,z)  =  L  +  5L« 
The  phase  difference  between  successive  wavefronts  is  then  given  by 
(w/c)6L.  Since  the  wave  proceeds  from  one  surface  to  the  next  in  time  6t 
while  the  phase  at  any  fixed  position  changes  at  the  rate  w,  this  phase 

difference  must  be  equal  to  w6t.  Further,  if  6   is  the  distance  between 

'     sn 

the  adjacent  surfaces  and  if  v  is  the  wave  velocity,  then  v5t  =  6   . 

sn 
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We   then  relate   these   expressions: 


(w/c)6L  =  loSt  ~(w6      )/v 

sn 


or  6L/6   =  c/v  ~r\    ,  index  of  refraction  of  the  medium  (2«8) 


Note  that  expression  2.8  is  the  directional  derivative  of  a  scalar  L 

in  the  direction  of  the  normal  to  the  wavefront.  By  definition  of  the 
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directional  derivative  of  a  scalar   function  4>, 


dc^/ds   =  T«V<> 

where  T  is  the  unit  vector  indicating  direction  in  which  the  derivative 
is  taken.  By  definition  of  the  vector  dot  product,  with  0  the  angle 
between  T  and  V<J)# 

d<j)/ds  =  T'V<J>  =  |T||v<j>|  cosO  -  |v<J>|cos0 

For  our  case  the  direction  of  interest  in  normal  to  the  front,  and  0=  0, 
Hence: 

6L/6   =  IVLl  -  c/v-  n  (2.9) 


llo 


Considering  the  curvature  of  rays  in  an  inhora  gei   us  medium  we  can 

demonstrate  that  the  rays  are  rectilinear  in  a  homogeneous  medium.  Let 
S  be  a  unit  vector  in  the  direction  of  a  ray  at  the  point  of  interest. 
This  vector  is  normal  to  the  wavefront  and  has  the  direction  of  VL.  Then, 

|vl|  =  p  or   |VL|  /n  =  1 

Then  S  =  VL  /n  (2,10) 

Let  N  be  a  unit  vector  in  the  direction  of  the  radius  of  curvature  of  the 

ray  at  the  some  point  of  interest  and  p  the  radius  of  curvature  of  the 
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ray  at  the  point.  The  vector  curvature  is  defined'  as  N/p  «  This  vector 

curvature  is  r.lso  given  by  dS/ds,  where  s  is  a  distance  measured  along 

the  ray.  By  definition  and  vector  identity  the  following  expansion  can 

be  made : 

dS/ds  -  (S«V)S  =  -  S  x  (V  x  S)  «  N/p  (2.11) 

This  expression  can  be  manipulated  to  show  a  simple  relationship  between 
curvature  of  the  rays  and  index  of  refraction  of  the  medium.  Details  arc 
left  to  appendix  1.  The  result  of  interest  is: 

1/p  -  N.V(ln  n)  (2.12) 


11 


Thus  for  a  homogeneous  medium  where  the  index  of  refraction  is  constant, 

or  independent  of  position,  ve  see  that  the  right  hand  side  of  1.12 
is  zero,  indicating  that  the  radius  of  curvature  of  the  ray  is  infinite, 
indicating  that  the  rays  are,  in  fact,  rectilinear.  It  follows  that 
VXS  r:  0}  which  is  a  sufficient  condition  for  the  existence  of  a  family 
of  surfaces  orthogonal  to  a  field  of  vectors  S. 

The  relationship  of  the  amplitude  of  fields  on  successive  \-7avefronts 
can  he  determined  from  the  premise  of  geometrical  optics  that  rays  are 

lines  of  flow  energy  and  that  no  power  flows  across  the  sides  of  a  tuhe 
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of  rays.  These  results  were  verified  in  the  prececding  paragraph.  Silver 

considers  two  surfaces  of  waves  L.  and  L  and  a  tuhe  of  rays  that  cuts 
out  elements  of  area  dA  and  dA  on  the  respective  surfaces.  By  con- 
servation of  energy  requirements,  the  flow  across  any  section  nortaal  to 
the  tube  will  be  constant.  In  terms  of  the  Poynting  vector  defined  by 


S  =  %(EXH  ) 


we  can  formulate  the  flow  constraint  by 


Sl  dAl  =  S2  dA2 


Since  in  free  space  the  Poynting  vector  is  proportional  to  the  magnitude 
squared  of  the  field,  between  successive  wavefronts  the  following  relation 
will  hold: 

3L8 


|E. |2dA.  -   |eJ2  dA„  (2.13) 


•xi  -"! 


Considering  the  ray  through  poiiit  A  on  the  surface  L  coincident  with 
the  z~axis  in  figure  5,  let  the  xzsy;j  planes  coincide  with  the  principle 
planes  of  L.  at  A,  A  ray  through  an  adjacent  point  B  lying  in  L.  and  the 
x-z  plane  will  intersect  the  ray  through  A  at  the  point  0  ,  a  distance 
R.  away.  This  distance  R.  thus  is  one  of  the  principle  radii  of  L.  at  A. 
A  similar  analysis  reveals  the  definition  of  the  other  radii  R„  based 
upon  intersection  of  rays  passing  through  adjacent  points  in  the  x-y 
plane.  The  point  A'  on  L_  lies  on  the  ray  through  point  A.  At  A1  local 
coordinates  x',  y1  are  constructed.  Since  the  fronts  are  spaced  apart 
by  the  distance  P,  the  two  principle  radii  of  curvature  of  the  wavefront 
L„  at  the  point  A'  are  R  +  P  and  R  4-  P. 

Considering  the  elements  of  area  dA.  and  dA  around  the  points  A,  A' 
bound  by  the  curve  Cf  C'  on  surfaces  L.  and  L_  respectively,  we  can 
write: 


dA  -  /(xdy-ydx) 

C  (2.14) 


dA.  -  //x'dy'-y'dx') 
2   c!r 


From  inspection  of  figure  5,  and  by  the  law  of  sines  we  can  relat; 
various  distances: 
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x1  -  (R.  +  P)x 


R„ 


(2.15) 


By  relating  dx,  dxf,  dy,  dy1  and  substituting  into  1.14  we  obtain  a 
relationship  between  the  area  elements: 


dA2  =  dAi  (Rx  +  ?) (R2  +  P) 


R1R2 


Substitution  of  2«16  into  2*13  we  obtain: 


(2.16) 


|E.|2dA.  -  |E9|2dA.  J&,  +  P)(R?  +  P) 


R1R2 


For  the  case  of  a  wavefront  diverging  from  a  region  of  sources  this 
expression  relates  the  amplitude  at  the  wavefront  farther  from  the 
source: 


3c 


E2  "  El 


_£*  R? 


(2-17) 


From  the  prcceediiig  development  it  would  appear  that  geometrical 
optics  is  based  on  the  local  behavior  of  the  wave front  as  a  plane 
wave.  Indeed,  this  is  the  basic  tenet  of  the  theory.  The  plane  wave 
solution  to  the  wave  equation  may  be  used  to  show  that  in  the  limit  of 
decreasing  wavelength  the  geometrical  optics  solution  satisfies  the 
Maxwell  fields  equations.  The  vector  solutions  appropriate  to  describe- 
linearly  polarized  electromagnetic  fields  that  behave  locally  like 
plane  waves  have  the  form: 

E  -  A(x,y,z)  eKo)t-koL(x,y,Z)) 

(2.18) 
H  -  B(x,y,z)  eKo>t-koL<x,y,Z)) 


The  amplitude  vectors  may  be  complex,  but  phases  must  be  independent  of 
position  to  fit  the  plane  wave  solutions.  Substituting  these  forms  of 
solution  into  the  homogeneous  forms  of  Maxwell's  equations  given  by 


V  X  E  +  jwuH  a  0 

(2.19) 
-jweE  +V  X  H  -  0 
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end  expanding  the  resulting  expressions t    these  may  be  solved  for  the 
vector  amplitudes  A  and  B  : 


A  ~^}±   (  VL  x  B)  +  1  (  V  x  B) 

o.»e  jus 


B  ^ k  (  VL  x  A)  —  _1_(  V  x  A) 


(2.20) 


where  k  is  the  free  space  propagation  constant,  k  -  w/ut   , 
o  o      o  o 

Equations  2„20  may  be  solved  to  yield  an  expression  in  terms  of  A  alone. 
Using  the  definition  of  k  and  the  index  of  refraction 


n  »  c/v  »   (eu/e  uq) 2 


we  obtain  the  following  expression: 


A  =  ^  (  VL  x  VL  x  A)  •!•  _1_(  VL  x  V  x  A  +  V  x  VL  x  A)  + 

n2"  Jkon? 


^  V  X  7  X  A)  (2.21) 


Expanding  2.21  with  the  vector  identity 

VL  x  VL  x  A  ^   L(VL'A)  -  A(VL'VL) 

we  obtain  for  A  : 
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A  -  ^l^  {  VL(A'VL)  -  A|VL|2}  +  _1_  (  VL  x  V  x  A  -:•  V  x  VL  x  A) 
r\2  jk  n2 


+  X   V  x  V  x  A  (2*22) 

o 


A  similar  expression  results  for  B,  replacing  A  with  B  in  2.22.  If 

L  and  the  first  and  second  derivatives  of  A  and  B  are  finite,  the  last 

two  terras  of  2.22  are  of  the  order  1/k  and  1/k   compared  to  the  first 

o        o 

term.  As  the  wavelength  approaches  zero,  k  approaches  infinity,  and 
we  are  left  with  the  following  expressions: 


A  ■  Zi-  {VL(A«VL)  -  a|vl|2  } 

(2.23) 


n 


B  =  -1_  {VL(B«VL)  -  B|VL|2  } 


For  these  expressions  to  yield  identity  of  amplitude  of  the.  vectors 
A  and  B  we  see  that  the  following  must  be  true  for  the  zero  wavelength 
case: 

A«VL  '■-   0 

B'VL  =  0 

These  two  conditions  note  that  A  and  B  are  transverse  to  the  gradient  of 
L,  or  that  they  both  lie  in  a  plane  transverse  to  the  direction  of  propa- 
gation. Furthermore,  since  2.20  expresses  B  in  terms  of  A  we  have: 
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B  »  /^7IT{(7L  x  A)  ~  JLJ7  x  A   }  (2,24) 

n        jkon 

We  sec  Chat  in  the  zero  Wavelength  case  the  second  term  vanishes 
and  we  are  left  with  an  expression  that  implies  vector  amplitude  B  is 
perpendicular  to  A  as  well  as  transverse  to  the  direction  of  propagation, 
since  VL/n  is  a  unit  vector  in  the  direction  of  propagation. 

Thus  the  field  vectors  of  geometrical  optics  possess  the  properties 
of  plane  waves  in  what  might  be  called  the  far-zone  fields,  indicating 
great  distances  from  the  source  with  respect  to  wavelength,  The  condition 
for  neglecting  terras  of  the  order  1/k  and  1/k  for  short  wavelengths 
require  that  all  associated  derivatives  are  finite.  In  the  neighborhood 
of  geometrical  focal  points  or  caustic  surfaces,  the  function  L  varies 
rapidly,  and  geometrical  optics  predicts  infinite  amplitudes  ,  which 
must  be  rejected  on  non-physical  grounds.  At  the  region  of  geometrical 
shadow  the  approximations  begin  to  predict  erroneous  results  due  to 
the  rapid  variation  of  amplitude  over  distances  small  compared  to 
wavelength.  In  these  regions  the  simple  methods  do  not  work,  and  dif- 
fraction and  scattering  theory  must  be  used, 
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Further  work    has  determined  that  the  Huygcns-Frcsncl  relation  in 

the  zero-wavelength  case  approaches  the  point-to-point  amplitude  relation- 
ship predicted  by  the  geometrical  optics  concept  of  energy  flow  in  tubes 
of  rays.  Agreement  holds  in  the  far-zone  approximation,   at  distances 
greater  then  wavelength,  and  reveals  the  constraint  on  amplitude 


W 


A's  normal  derivative  that,  besides  being  finite  everywhere, 


JL  J>A  «  2tt 
A  3n"   '  \ 


This  is  satisfied  in  the  limit  of  zero  wavelength  providing  the  frac- 
tional change  in  amplitude  over  distance  equal  to  that  of  a  wavelength 
is  small  compared  with  unity. 


Since  the  basic  tenet  of  geometrical  optics  requires  that  the  wave- 
front  behave  locally  like  an  infinite  plane  wave,  scattering  phenomenon 
may  be  interpreted  geometrically,  and  appropriate  laws  of  reflection 

derived  from  the  field  equations  and  boundary  conditions.  That  these 
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results  are  valid  can  be  developed  from  other  principles  of  optics 

including  Fermat's  principle  and  the  law  of  the  optical  path,  demonstrat- 
ing the  properties  of  optical  rays  and  the  reflection  phenomenon,  or 
Snell's  laws. 

Me.   consider  the  two-dimensional  analysis  of  an  infinite  plane 
wave  incident  upon  an  infinitely  conducting  (  that  is  to  say  perfectly 
reflecting)  sheet.  The  geometry  of  the  problem  is  described  by  figure  6. 
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6^s>  CO 


"Ti^iUV-e  (o 


The  incident  fields  will  have  the  form: 


.   ,    -ik  z  cos0.  +  jk  x  sin0. 
-  i  h.  e  *'  o      i   J  o      i 


y  i 


(2.25) 


_   „/.      ,     „^„~      i     •      u      •    ^   \      ~jk   z   cog0.    +  jk  x  sinO. 
=   £(i     h.cosQ.   +  i     h  sm0.)    c  J    o  i    •    J    o  i 


7-     i< 


v;here  5    is   the  wave   impedance   of    free   space,    equal   to    (^    /e   ) 

o  o 

The  reflected  fields  will  have  the  following  form; 


tI  _  .   ,    +jk  z  cose  +  ik  x  sinO 
H  =  1  h  e  J  o      r     o      r 
r    y  r 


(2.26) 


Er  =    rf"1!  hrcos0r   +  i_  h„sinO„)e^koZ   COr:°r  *  jV  sinGr 


r         z     r         r 


3b 


Since  at  Z  =  0  Eor  all  ,  lmd  ,   the  bQundary  condiWon  on  thfi  tansentiai 

component  of  electr-in   f-f«u     ^ 

exectiic  field  requires  that   (E  +  E  )      te  0 

Then,  equating  reflected  and  incident  f-ui  i       '   ^^ 

a  incident  fields  appropriately  we  obtain: 


e  =  0 

i    r 


h.  -  h 
i    r 


c  »  e 
i    r 


Thus  we  note  th-it-  ft,-  *  -n 

that  the  Mta,  boundary  conditions  exist  at  the  inter- 
face: 


tang 

(2.27) 


x  nor    v  r' 


r  nor 


These  expressions  may  be  given  in  terms  of  a  Vttfnr   ,   , 

s  oi  a  vector  relationship: 


n  x  (  Ej  +  E  )  =  0 

„  .*      „  (2.28) 

1      r 

We  can  examine  this  vector  r.i.M   i  • 

vector  relationship  to  determine  an  .„ 

j     "cj.erra5.ne  an  expression  for 

-  «-*  faceted  fron  .  sor.ee  in  terns  of  the  ^  ^  ^ 
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With  the  vector  form  of  the  boundary  conditions  we  form  the  cross  product 
of  the  normal  n  and  the  first  expression  of  2.28e  We  then  multiply  the 
second  of  2.28  with  the  normal  vector  ri  and  add  the  two  expressions. 
This  yields; 


n  x  E  x  n  -r-  (n*E  )n  =  -n  x  E.  x  n  +  (n«E.)n 
r  r  i  i 


Expanding  by  the  vector  identity 


A  x  B  x  C  «=  (A«C)B  -  (A»B)C 


we  obtain  the  expression 


E  «  -E.  +  2(n«E.)n  (2.29) 

r     1       i 


This  may  be  interpreted  to  express  quantitatively  the  fact  that  the 
reflected  field  experiences  reversal  of  its  sense  of  tangential 
components  in  order  to  provide,  a  zero  tangential  field  at  the  surface 
of  the  reflector. 


Thus  we  have  demonstrated  that  the  approximate  methods  of  geometrical 
optics  involving  point-to-point  transformations  along  rectilinear  rays 
between  successive  wavefronts  or  surfaces  of  constant  phase  can  provide 
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adequate  phase  cud  amplitude  information  providing  certain  constraints 
are  observed.  Furthermore,  that:  thcae  fields  satisfy  Maxwell's  fields 
equations  and  the  wave  equation  within  the  limits  of  certain  constraints 
has  been:  demonstrated  by  a  zero-wavelength  limit.  We  summarize  the  results 
of  geometrical  optics  in  which  we  arc-,  interested  as  relevant  to  this  works 

1.  Energy  passing  through  consecutive  cross- 
sections  of  a  tube  of  neighboring  rays  is 
constant, 

2.  Rays  in  a  homogeneous  medium  are  rectilinear. 

3.  At  a  surface  of  reflection,  an  incident 
wavefront  behaves  locally  like  a  plane 
wave,  and  the  reflection  of  rays  may  be 
geometrically  interpreted  in  accordance 
with  specular  reflection  laws. 

The  following  constraints  must  be  observed  for  the  geometrical  optics 
approximations  to  holds 

1.  All  lengths  of  interest,  radii  of  curvature, 
for  example,  must  be  large  compared  to  a 
wavelength. 


33 


2.  Spacing  between  neighboring  rays  must 
change  very  little  over  distances  on  the 
order  of  wavelength. 

3,  The  fractional  change  in  amplitude  of 
field  quantities  over  distances  on  the 
order  of  wavelength  must  be  small  com- 
pared to  unity. 

A.  Geometrical  optics  cannot  be  used  in 
the  region  of  geometrical  shadow. 

5.  Geometrical  optics  fails  at  focal  points 
or  caustic  surfaces,  predicting  non- 
physical  results. 


ho 


2.4   It  remains  to  utilize,  the  previously  developed  expressions  to  describe 
the  phase,  amplitude  and  direction  of  the  electric  field  over  the  ref- 
erence sphere.  The  choice  of  radius  of  the  reference  sphere  is  discussed, 
as  is  the  technique  of  ray  tracing  to  be  used.  Expressions  for  phase 
and  amplitude  on  the  reference  sphere  are  developed,  relating  thase 
quantities  to  those  incident  upon  the  reflector  surface.  Finally,  the 
distribution  over  the  reference  sphere  is  related  to  a  coordinate  system 
with  coincident  origin. 

For  the  geometry  of  figure  7,  we  consider  determination  of  the  fields 
produced  by  some  source,  directed  towards  the  reflector, that  will  produce 
illumination  over  the  aperture  in  the  form  of  waves  propagating  in  the 
-  z   direction.  The  fields  propagating  in  this  direction  will  be  of  the 
form: 


E,  H  «  eJ  o 


or  in  spherical  coordinates: 


„  TT     ik  R  cosG 
E,  H  «  e  o 


Since  we  desire  to  specify  these  fields  over  a  reference  sphere  of  radius 
R'  with  origin  at  z  =  cR,  where  c  is  some  constant  fixing  the  origin, 
the  phase  expression  must  be  corrected  to  account  for  the  additional 
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+z>  referee 
sphere 


V* 


rfliM  1 


~XU  1 


incident  ra^ 


path  length  of  travel  to  the  reference  sphere.  Inspection  of  figure  7 
Indicates  that  the  phase  variation  has  two  cases  of  interest,  one  for 
angles  o  >  0  auc^  'ne  other  for  the  case  of  6  <  G  ,  »  where  0  is  the 
angle  rt  which  the  reference  sphere  intersects  the  spherical  reflector. 
The  variation  in  field  expressions  is  introduced  as: 


r>   it        jk   (Rcos  -5-  L, ) 
G  >  ec    E,H       eJ  o  V 


r  tr        jk   (Rcos  -  L„,) 
0  <  e     E,H   «    eJ  o  2/ 

c 


(2,30) 


Inspection  of  the  exponent  suggests  that  if  appropriate  values  of  other 
parameters  could  be  found  so  as  to  keep  the  difference  to  a  minimum  that 
synthesis  might  be  enhanced,  since  the  smaller  exponent  represents  a 
slower  varying  sinusoid.  Furthermore,  amplitude  considerations  in  the 
geometrical  optics  approximation,  namely  the  desire  to  avoid  the  region 
of  the  caustic  surface,  suggest  that  a  vise  choice  vould  be  to  consider 
the  reference  sphere  intersecting  the  reflector  surface  at  the  latter' s 
extreme  edge.  Thus  the  radius  of  the  reference  sphere  will  be  constrained 
by  the  radius  of  the  reflector  under  investigation,  R,  the  position  of 
the  origin  of  the  reference  sphere  z  =  cR,  and  the  maximum  angle  sub- 
tended by  the  reflector.  By  laws  of  trigonometry  we  observe: 


R'2  =  R2  +  c2R2  -  2cR2cos6  „  ,n 

m  (2.31) 
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The  ray  tracing  technique  proceeds  from  a  ray  incident  upon  the 
reflector  at  a  particular  angle  of  interest*  The  ray  reflects  in  ac- 
cordance with  the  previously  developed  laws  of  reflection.  We  then 
assume  that  if  a  hypothetical  source  could  produce  the  appropriate 
fields  that  would  give  the  same  illumination  as  we  have  upon  the  aper- 
ture, these  rays  would  remain  as  we  have  determined,  but  the  sense  of 
propagation  of  energy  along  the  rays  would  be  reversed.  Vie  then  remove 
the  reflector,  and  trace  the  rays  until  they  intersect  the  reference 
sphere.  Within  the  limits  of  geometrical  optics,  v;e  can  estimate  the 
amplitude  and  plicae  at  each  point  of  intersection,  based  upon  the 
extra  path  length  travelled  by  thu  waves,  along  the  rays,  to  the.  refer- 
ence sphere. 

It  is  necessary  to  calculate  the  distance  L.  We  proceed  with  an 
equation  determining  the  distance  between  points  on  the  reflector  and 
points  on  the  reference  sphere,  as  measured  along  the  ray  of  interest. 
Referring  to  figure  8: 


L2  -   (x  -  x  )2  +  (y  -  y  )2  +  (z  -z  )2     (2.32) 

a    m      Ja       J\n  am 


where 


x  -  R'  sinO'  coscf)'    x  =  R  siuG  cos<J> 
a  m 


y  -  R'  sinO'  sin*'    y  -  R  sinO  sin* 


z     «=  R'  co^O'  +  cR     z   =  R  co. 
a  m 
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Noting  from  symmetry  that  <}>  s  <*>',  and  substituting,  expanding  terms 
we  obi r.in: 


(  .         L  =    (R'sinO'  -  RsinO)2  +   (R'cocO1  4  cR  -  RcosG  )2 


To  proceed,  vze  must  relate  0  and  0'  .  Vievjing  in  the  x-z  plane,  the 
equation  of  ray  is  given  by: 


x  -  (z  -  z_)  tan2 
y  -  0 


°  (2.3/») 


From  the  equation  of  a  spherical  surface  : 

R'2  =  x2  +   (z  -  cR)2  +  y2  (2.35) 

We  can  solve  2.34  and  2.35   in  the  plane  y  =  0  to  yield  an  expression 
representing  the  intersection  of  a  ray  and  tho  surface  of  the  reference 
sphere.  Solving  for  z  ,  intercept  of  the  z-axis  we  obtain: 

R    E    z 
sin7^~2e)  ""sinO  (2.36) 

Substituting  2.34  into  2.35,  and  expanding  in  a  straightf orward  fashion 
we  obtain,  with  substitution  of  2.36: 


z2  -  2z(  Rsin0sin20  +  cRcos226  )  4  R2sin26  -  cos220(R*2  -  c2R2)  -  0 

Details  of  the  manipulation  are  included  in  appendix  1,  Normalizing  all 
distances  to  Rt  radius  of  the  reflector  under  investigation,  we  obtain: 


z  2»  2z  (  sin6sin29  +  c  cos220  )  +  sin20~  cos220(R'2  -  c2)  -  0 
n     n  n 

(2.37) 
where 

R1  =   R'  z  =   z 

n   __  n   __ 


Solving  the  quadratic  2.37  by  usual  methods,  we  obtain  for  z      : 


z  =  (  sin0sin29  +  c  cos2 20)  + 
n 


(sinOsin20  +  c  cos220)  -  sin20f  cos220(R*2  -  c2) 

n 

(2.38) 


Noting  from  the  geometry  of  figure  8: 


R'cos    '   -  z'   «  r-  cR 


or 

cos0'      -      (z  -   cR)        -        (z     -  c)  (2.39) 

n 


4k 


Substituting  equation.  2.38  into  2.39  the  foil   '    relationship  betv. 
angles  8  and  G'  is  obtained: 


cone'  -  1   (  sin0sin20  +  c  cos229  -  c 

R» 

n 


+   (sin6sin29  +  c  cos229  )  -  sin2e  +  cos229(R'2  -  c2)      ) 

n 


(2.40) 


Equation  2.40  states  that  there  is  a  relationship  between  angles  such 
that  for  every  angle  of  incidence  on  the  reflector  surface  there  exists 
an  angle  0*,  as  measured  from  the  origin  of  the  reference  sphere,  that 
describes  the  intersection  of  the  reflected  ray  traced  back  to  the  refer- 
ence sphere  and  the  surface  of  this  sphere.  Thus  the  angles  may  be 
related  in  terms  of  the  parameters  of  the  reflector  under  investigation, 
R  and  9   ,  and  the  center  of  the  reference  sphere  implied  by  the  value 
of  c.  Once  the  angles  associated  with  a  particular  incident  ray  have  been 
related,  they  may  be  substituted  into  2.33  to  yield  the  distance  L  as 
a  function  of  the  point  of  interest  on  the  reference  sphere. 

Equation  2.17  was  developed  to  express  the  amplitude  of  field  quanti- 
ties on  one  wavefront  with  respect  to  the  field  quantities  on  another 
front.  It  is  seen  that  the  amplitude  varies  with  the  two  principle  radii 
of  the  wavefront  and  the  distance  of  separation  between  fronts.  It  is 
necessary,  then,  to  develop  expressions  for  the  radii. 
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Once  again  we  consider  the  case  of  a  spherical  reflector  illuminated 
by  aa  incoming  plane  wave*  A  typical  reflected  wavefront  has  been 
constructed  in  figure  9 ,  where  the  froxit  represents  a  surface  of 
constant  phase  as  defined  by  the  total  path  length  measured  along 
each  ray,  L  +  L„  K  R  .  The  wavefront  is  formed  by  connecting  normals 
erected  on  each  ray  at  the  point  where  total  path  length  equals  the 
radius  of  the  reflector.  Figure  9  shows  the  front  as  viewed  in  the 
plane  y  n  0.  From  symmetry,  it  may  be  seen  that  the  surf nee  of  the 
front  can  be  generated  by  revolution  of  the  projection  in  figure  9 
around  the  axis. 

At  any  point  of  interest,  a  wavefront  may  be  passed  through  this 
point  by  choice  of  the  appropriate  path  length.  At  each  point  the 
wavefront  will  have  two  principle  radii  of  curvature.  Since  a  ray  "is 
normal  to  the  front  and  any  two  adjacent  rays  intersect  at  the  caustic 
surface,  forming  the  envelope  that  is  the  caustic  surface,  it  is  seen 
that  one  radius  of  interest  is  the  distance  from  the  front  at  the  point 
of  interest  to  the  point  on  the  caustic,  as  measured  along  the  ray.  The 
second  radius  of  curvature  may  be  determined  from  the  axial  symmetry  of 
the  system.  In  the  <{>  direction,  adjacent  rays  will  intersect  at  the  axis 
of  symmetry,  and  the  second  radius  of  curvature  will  be  the  distance 
from  the  front  at  the  point  of  interest  to  the  axis,  measured  along  the 
ray. 
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Rays  incident  upon  the  reflector  at  an  angle  0  are  reflected  at  this 
angle,  pass  tangent  to  the  caustic  surface  and  cross  the  axis.  The 
coordinates  of  the  point  to  which  tangency  eccurs  are  given  by  equations 
2,6  and  2.7,,  the  parametric  equations  of  the  caustic.  From  figure  9  we 
see: 


R.  =  MU  =  AT~~~  z  yT-'jr—Q     „   x  JT 
1  m    c        n    c 


(2.A1) 


R2  =  MW  - 


z  -  z  Jz~T~{k     -  W 

no        m 


(2.42) 


where 


z  ~ 
ra 


x  = 

m 


Rcose 
RsinS 


z  =  R 


sin 


sin2e 


Substituting  into  equations  2.41  and  2.42,  and  normalizing  all  quantities 
to  the  radius  R  we  obtain: 


R 


In 


COSQ 

2 


(2.43) 


R 


2n 


(2.44) 


2  cose 


Details  of  the  manipulation  involved  are  covered  in  appendix  1. 


So 


Although  the  wavefront  has  been  analyzed  by  considering  illumination  of 
the  reflector  by  an  incoming  wave,  it  is  the  fact  that  if  a  source  can 
be  provided  that  can  generate  such  a  wavefront  as  has  been  described,  the 
wavefront  will  be  reflected  from  the  spherical  surface  to  yield  an 
outgoing  plane  wave  over  the  reflector  aperture.  Figure  9  represents  a 

graphical  illustration  of  the  method  of  compensation  suggested  by 
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Ashmead  and  Pippard,  a  vavefront  whose  departure  from  sphericity  is 

just  that  needed  to  transform  into  a  plane  wave  upon  reflection. 

Since  we  wish  to  describe  the  fields  over  the  reference  sphere  as 

a  function  of  coordinates  with  origin  at  that  of  the  reference  sphere, 

wc  must:  introduce  a  transformation  of  the  field  vectors.  Initially, 

field  quantities  were  described  upon  incidence  in  terms  of  coordinates 

with  origin  at  that  of  the  reflector.  The  geometry  of  interest  is 

described  in  figures  10  and  11,  considering  both  cases  where  the  center 

of  curvature  of  the  reference  sphere  is  exterior  and  interior  to  the 

intersection  of  the  reflected  marginal  ray  (edge  of  reflector  being  the 

point  of  incidence)  and  the  axis.  Note  that  due  to  axial  symmetry  there 

is  no  modification  necessary  to  the  (}>  components  By  inspection: 


i  ,  »  i  cos  (G'  -  0)  +  i  sin  (0*  -  6) 
r    r  G 


iot  «=  -  i  sin  (01  -  6)  +  iQ  cos(Q'  -  0) 
or  o 


(2.45) 


Si 


</ 


•A —    i<\ci4on"t  taw 


TiQut-e  i-° 


W 


doxctul 


\ociacvH-  v-av 
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Solving  for  the  imprinted  components   in  terms  of   the  primed   system: 


lr  -   i^cos  (0f    -    0)    -   iG,sin(0*    -   0) 


ie   =   irtsi.n(0f    »   e)   +  iotcoc(0'    -   0) 


(2.46) 


From  the  preceoding  derivations  it  may  be  seen  that  equations 
2.17,  2c 30  and  2.29  provide,  sufficient  information  to  form  an  expression 
for  the  geometrical  optics  amplitude,  phase  and  direction  of  the  electric 
field  distribution  over  the  reference  sphere,  based  upon  a  desired 
reflected  field  E   .  Quantities  of  interest,  including  the  radii  of  curv- 
ature of  the  vavefront,  distance  from  the  point  of  the  reference  sphere 
under  consideration  to  the  corresponding  point  on  the  reflector,  and  the 
relationship  between  angles  are  given  by  2.43,  2.44,  2.33  and  2.40 
respectively.  Finally,  transformation  of  the  vector  components  given  by 
2.46  relates  this  field  distribution  to  coordinates  at  the  center  of  the 
reference  sphere.  Thus  the  geometrical  optics  field  is  described  over  the 
reference  sphere  as  a  function  of  the  reflector  surface  under  investiga- 
tion and  the  particular  origin  of  interest  of  the  reference  sphere. 
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In  summary,  we  formulate  an  expression  for  the  general  ease  where 
a  desired  electric  vector  field  is  specified  over  the  aperture.  Choosing 
this  field  to  have  unity  amplitude,  we  describe  the  field  that  would 
occur  over  the  reference  spheres 


Amplitude!   E  c  =  E   ,. 

r  reference     reflector 

sphere 


Ri  Ro 
In  2n 

(R,   *  L  )  (R.  +  L  ) 
in    n   2n    n 


*2(2.17) 


Phase: 


r»i     /  r,  \  ik  R(  cosO  -  L  ) 

Phase (  E  ,      )   «  eJ  o  v         n' 
reference 

sphere 


(2.30) 


Direction;   E 


reference 
sphere 


-  E  +  2(  n  -E  )n 
r         r 


(2.29) 


Combining: 


E 


reference 

sphere 


......  R^R?\. 

(R   +  L  )(R   +  L  ) 
In    n  la         n 


h     jk  R(cos 
cJ  o 


V  {~E  -:•  2(n-E  )n) 
r       r 


(2.47) 


SH 


where: 


R,     cosG 
In  =  — - r— 


(2. A3) 


R 


2n 


1 

2cosG~ 


(2.44) 


R' 
n 


=  1  +  cz  -  2c  cos( 


(2.31) 


n 


(sinG  -  R'sine')   +  (R*  cosG'  +  c  -  cosO)  (2.33) 
n  n 


cose'  ~  -Jr—  (  sin0sin28  -I-  c  cos220  -  c 

R' 

n 

+   sinGsin2G  +  c  cos226)2-  sin29+  cos22G(R'2  -  c2)  )  (2.40) 


i  -  i  ,003(6'  -  6)  -  i.fsin(9'  -  G) 
r    r  v 


i  -  i  ,sin(e'  -  G)  •;-  i  ,cos(6'  --  G) 


(2.46) 


SS 


} 


Twu^e  2-2- 


E.=  k  'V  e 


J 


K»3 


aperture,  i\<\e  paction 


2.5   In  this  section,  we  apply  the  results  of  the  general  expression  2.47 
to  a  particular  case  of  interest.  Specifying  the  field  at  the  aperture 
shown  in  figure  12  to  be  a  uniform  plane  wave  linearly  polarized  in  the 
x-  direction  we  have: 


aperture 


jk  z 

1   cJ  o 

X 


Expanding  in  spherical  coordinates: 


E  =     eJk0Rcos3  sineCos<J>   +  i„cos9cos<{>   -   i.sin*)         (2.48) 

aperture  r  6  T<j)r/V/ 


By  applying  equation     2.29t    reversing   the  sense  of   tangential   components, 
we  can  determine   the   sense  of  direction  of   the   field   over   the   reference 
sphere: 
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direction(  E   -      )   ~  i  sinOcosA  -  i^cosOcor;A  +  i.sindi   (2.A9) 
reference       r  6  $ 

sphere 


Applying  equations  2.46  to  2.49  the  direction  of  the  field  vector 

cr.n  bo  expressed  in  terms  of  the  origin  at  that  of  the  reference  sphere: 


I  -  i  ,cos(e'  -  e)  -  i, since1  -  e) 


i0  -  ir,sin(G'  -  6)  +  i^cosCO'  -  0) 


(2.46) 


Substituting; 


direction(  E   _      )   =»  i  ,cos$  (  sinOcos(0'  -  0)  -  cos0sin(0'  -  0)) 
sphere 


-  i  »cos<{)(  sin0sin(0'  -  0)  +  cos0cos(0'  -  0)) 


+  ixsin<f>  (2.50) 


Collecting  terms  and  applying  trigonometric  identities  for  the  difference 
of  angles  we  obtain: 


direction  (  E   ,      )   =  i  fcos*sin(20  -  0C)  -  i  ,cos<>cos(2G  -  0') 
reference       r    Y  0 

sphere 

+  i,sinc(>  (2.51) 

S~7 


Thus  ft-  previously  developed  relationshlps  2.17  and  2<30>  .^..^ 
for  amplitude,  phase,  and  direction  of  the  electric  field  over  the  ref- 
erence sphere  can  be  combine,  to  yield  the  vector  field  solely  as  a 
function  of  the  particular  reflector  under  investigation  and  the  choice 
of  .rigia  of  the  reference  sphere  along  the  axis. 

While  amplification  could  be  perf.mee  on  the  rcsluting  expression, 
the  field  over  the  reference  sphere  is,  in  itself,  tnl>.  ,  „„  t# 
eventually  describe  the  fields  over  a  snaller  sphere  surrounding  s-e 
source  region,  concentric  with  the  reference  sphere.  Since  the  sta.lifica- 

tion  would  not  enhance  this  procedure  4+   „*n    ,    i 

prsceaure,  it  will  net  be  pursued.  Rather, 

the  resulting  expression 


"reference 
sphere 


In   2n 


(R 


ln         n        2n  n 


H         JkR(   cose   »  L   ) 

CO  n 


(   i    ,cos$siii(2e   -   6')    -   i 


,cos<^cos(2e  -  9*)  +  i  sind.  ) 


(2.52) 


(where  R^,  R^f  ^  and  anglfi  relatlonnblps  arfi  deflned  by  2  ^ 

2.44,  2.33,  and  2.40  and  are  listed  for  convenience  in  section  2. A) 


lends  itself  readily  to  tecniques  of  machine  computation.  As 


the  following  pages  ccnta 
in  amplitude  and  phase  fo 


a  result, 


in  graphs  representing  the  field  component.- 
r  a  particular  choice  of  parameters.  Thr 


*re  calculated  for  normalized  amplitude  and  reflect  no  variat< 

S9 


asimuthal  i  directiaa.  Thus  the  components  E  /  and  E   indicated  in  the 

u       <J> 

graph  must  be  Multiplied  by  cos$  and  r>in§   respectively  to  represent  the 
true  nature  of  the  field  over  the  entire-  sphere.  Similarly,  the  E  ,  cem- 
peneat  must  be  multiplied  by  the  variation  cost*;  ,  The  computer  program 
developed  to  calculate  the  field  over  the  reference  sphere  is  included 
in  appendix  II. 
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CHAPTER  III 
ELECTROMAGNETIC  FIELD  EQUATIONS,  EXPANSION,  AND  SOLUTION 

3,1   In  the  preceeding  chapter  expressions  for  the  electric  field  over  the 
reference  sphere  x^ere  determined.  In  this  chapter  those  fields  will  be 
matched  at  the  reference  sphere  to  an  expansion  of  spherical  harmonics 
following  from  solution  to  the  vector  wave  equation  in  spherical  coordi- 
nates. These  harmonics  will  be  evaluated  at  a  particular  radius  within 
that  of  the  reference  sphere  to  yield  the  electric  field  over  the  surface 
of  that  smaller  sphere  enclosing  some  source  region.  These  fields  will  be 
compared  to  the  solution  of  Ricardi  for  the  same  region.  Finally,  these 
results  and  the  comparison  of  solutions  will  be  interpreted. 
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It  may  be  shown  that  solution  to  the  scalar  Helmholtz  equation 


(  ?2+  k2  )y  =  0  (3.1) 


for  spherical  coordinates  is  given  by  separation  of  variables  to  be; 


V    -  b.  (kr)  Pm(cose)  {C?9Tnt>  (3.2) 

mn     n      n        sinm<j> 


where,  b,  (kr)  is  the  spherical  Bessel  function  and  P  (cos9)  is  the 
n  r  n 

associated  Legcndre  function  of  the  first  kind.  The  choice  of  cosing 
or  sinm<£  depends  upon  the  variation  of  the  fields  at  the  boundary  of 
interest. 


bZ 


The  fields  under  investigation  in  chapter  2  were  shown  to  be  vector 
fields,  and  the  solution  to  the  vector  Helmholtz  equation  must  be 
considered  as: 

(V2  +  k2)A  ■  0  (3.3) 

where  A  is  a  vector  field,  representing  either  the  magnetic  vector 

potential  or  its  dual,  the  electric  vector  potential.  It  has  been 

35 

shown  that  a  general  solution  may  be  formulated  by  considering  a  set 

of  fields,  both  magnetic  and  electric,  which  are  transverse  to  the 

radial  vector,  i  .  In  this  case  the  solutions  to  3.3  are  found  to  be: 

'  r 


A  -  Ya  r  i 
r 


F  -  Yf  r  i 
r 


(3. A) 


where  ¥  is  the  solution  to  the  scalar  Helmholtz  equation,  with  the  super 
script  indicating  that  is  solution  to  the  vector  Helmholtz  equation 
for  the  magnetic  potential  A  or  the  electric  potential  F. 


G3 


The  most  general  solution  for  fields  in  a  source-free  homogeneous  region 
in  terms  of  vector  potentials  is  found  to  be   a  linear  combination  of 
fields  with  the  behavior  of  equations  3<4,  where: 


E  -  -v  x  F  -  juu  A  +  JL_  v  (V?  A) 


H=  y  x  A  -  jaeF  +  _1_  V(V.F) 


(3.5) 


Since  the  solutions  for  A  and  F  involve  the  product  of  r  and  V ,  it  has 

been  found  convenient  to  adopt  an  alternative  definition  of  the  spherical 
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Eessel  function   : 


H  (kr)  -  kr  h  (kr) 
n  n 

The  superscript  H  indicates  we  have  chosen  the  spherical  Hankel  function 
of  the  second  kind,  whose  behavior  indicates  propagation  outward  from  the 
origin,  as  applicable  to  this  problem.  Anticipating  further,  since  the 
fields  we  wish  to  expand  are  of  the  form  of  equation  2,52,  where  the  var- 
iation with  (j,  is  given  by  sin<j>  and  cos<+>  ,  ve   choose  m  =  1  .  Thus: 


A<2>  1  f    sin*    V 

A     -      z   a     H    (kr)      Pi(cos8)J   S      6   { 

r  n     n  n  {  cos  v  } 


F     -     Z  b     H    (kr)      l^Ccose)1,   S1Il(f    ( 
r  nn  n  {   cos*    } 

n 


W 


(3.6) 


Expanding  equation  3.5  in  component  form  ve  obtain: 


E    =      l     ,  a2      -:-  k2  :  A 

r      _  (  _  )    r 

jwe        9r2 


,  3F        ,  ,  32A 

rsiii9     3<J>  jioer       3r3-0 


_  .       3F  ,  .  32A 

E,   »       1         r       +         1  r 


<> 


r      86  jwcrein.e        3r3$ 


a2 

H     -       1     ,  ._ +  k2       F  ,0   ,. 

r        — -  (  —  )     r  (3.7) 

jwy        3r2 


i  3A        ,      ,  32F 

H„  =       1  r     -5-     1  r 


0 


rsinS     3<}>  jcour     3r3i 


,      3A  .  32F 

r     3G  jwursinS        3r8<> 


Referring  to   the   component   EQ   in  3.7,   we  note   that   <j>  variation  for 

o 

3F 

A  and  J_r  is  the  sane.  Hence  ve  choose  the  appropriate  variation  in 

equation  3,6.  We  note  further  that  the  function  P  (cos0)  is  identically 
zero  for  n  «  0.  The  series  need  begin  for  n  =  1,  and  may  be  written: 


GS 


A  (**  x 

A       -  I     a       H   (kr)    P    (cos6)    cosi 
r  R       n  n 

n-1 


F       »   E     b        li    (kr)    P    (cos9)    sin* 
r  ,     n       n  n 


(3.8) 


Substituting  into   the  appropriate   components   of   3.7   the  electric   field 
can  be  formulated    : 


E       "       A_'r    H   (kr)     P    (cose)    cos(f)     i 

jwer2 


,  Aft)  1 

b      H    (kr)      P    (core)    ecscj,   1 
n     n  n  r      6 


rs3.nt 


„   _.                                 ,1                            ,    <JH    (kr)        dP1(cose) 
(3.9)  +  »  a     cos<t<  _JL__      ^-IL-, -1       *■ 


jwer  dr  d6 


1  ,       TT    ,,     N       .     ,       dP    (cose) 

r  d6 


a  (?-i 


,    A   ,,1,        ftN        dH    (kr)      . 
a     sind)   P    (cose)  nv  i 

n  n  —_-«—,         ,*, 

dr 


6<b 


For  convenience  we  define: 


d  II  (kr) 
n 


dr 


A  (21 

k  __£  Vkr) 
d(kr) 


k  H  (kr) 
n 


and  since  vre  are  concerned  with  the  case  of  propagation  in  free  space 


k  =  k  =w/c  u 
o     o  o 


C  -  Awe 
o  o 


Substituting  the  relationships  from  above,  multiplying  through  by  j / j  =  1, 
and  including  the  summation  we  form  equation  3.10  from  3.9  : 


(3.10) 


E  - 


n=l 


4-a     n(n-fl)  „    ..     .    -1,         _N 

n  II    (kr)    P    (cos 6)    cos 6    i 

~— — — — — —  an  r 

jwer 


+  j^a     H    (kr) 
n     n 


.     .   F    (cosG)  .     .    _,!,        ..     . 

-ccc<J>  _d   jij_ _  ._/_  xr,  +  S3-T^   P    (cosG)    l 

d   8  e        sine     n  < 


+       b     II   (kr) 


n     n 


-1^ 


x       n1/  o\     •        i  ^     i      i    P    (ccsO). 

■cosjb      r    (cosG)    x.   +      sin  A   d      n  i, 

sine"      ,l  8  T? * 


67 


Introducing  tangential  spherical  harmonics  as  defined  by  Morse  and 

„  ^  i38 
Fcshback  : 


BC      -  1 

mn 


/n(n+l) 


.    dP'  (cos 8)    ,  lTtn,        -v     . 

C0G:,V *k_ ^L  ifi  ~  mslnmftP    (cosG)    i 


(3.11) 


C°     -  1 


irm 


^n~(n!:l) 


d0  sine 


±     T>™/  n\     1  -f         j.     <*p     (cog  6)  . 

mcosEKfr  P  (cosO)  ip  -  sinm<J>   n    __  i 
.sinO  dO 


Substituting  3,11  into  3*10  we  can  form  an  expression  for  the  tangential 

fields.  By  exploiting  certain  orthogonality  properties  of  the  tangential 

harmonics,  ve  can  derive  the  constants  a  and  b   in  3.10  from  the  fields 
*  n      n 

specified  at  the  boundary  of  interest,  In  our  problem,  it  v;ill  be  the 

tangential  fields  on  the  reference  sphere  that  will  be  specified.  Note. 

that  from  equation  3.10  the  two  constants  a   and  b   are  sufficient  to 

n      n 

describe  the  field  vector.  It  must  also  be  noted  that  development  of  all 
expressions  in  this  chapter  has  proceeded  for  the  case  where  the  origin 
of  the  sphere  for  which  boundary  conditions  were  given  coincided  vzith 
the  origin  of  the  coordinate  system  for  the  spherical  harmonics.  It  new 
becomes  necessary  to  stipulate  that  the  harmonics  must  be  expanded  about 
an  origin  at  the  center  of  the  reference  sphere.  To  adapt,  these  expression 
we  need  only  to  substitute  0'  for  0,  recognizing  that  the  angle  is  de- 
fined as  in  chapter  2„  and  for  the  radius  variable  r  v/e  substitute  r', 
indicating  that  it  is  measured  from  the  prime  coordinate  origin. 
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From  the  preceeding  we  can  write  for  the  tangential  field  components: 

,       »    ^ r    At?-"»  a(z\ 

E            =     -  Z   /n(n+l)    lb     H    (kr)    C?       +  i^a     H    (kr)    B. 

tang  — -  Inn  In  n     n  In 

r         n=l  L 

(3.12) 


Coefficients  are  evaluated  by  forming  the  vector  dot  product  C,  »E    , 

la   tang* 

multiplying  each  side  of  3.12  by  sin9dGd.4>,  integratiiig  over  the 

interval  0-tt  end  0-2tt  on  $   and  0  respectively,  evaluating  the  fields  at 

the  reference  sphere,  radius  R' ,  The  following  properties  are  listed 
below: 


b!:    »B^  2TTn(n-!-l)  n=q 

/      /  sin0dGd<J>        -  2n+l 


lq      1ft 


n^q 


/  /   cJn*B?   sinGdGd<f>   =      0    for  all  n,q 


Applying  these  properties  to   equation  3.12  we  have: 


£9 


i(2n-!-l)    R'    //   E      «B®      sinOfde'd(J) 
a       ~     ___,_    „,,._ _  r_t     In 

C2ttH  (kR1)    (n(n+l))3/2 


-(2a+l)    R'    //   E      »C®     sine* dO1 d< 
b        =  rt      In 


(3.13) 


n 


2ttH  (kR')  (n(n-fl))3/2 


where  E   is  the  vector  representation  of  the  tangential  components  ©f 

the  electric  field  over  the  reference  sphere  at  radius  R1 «  To  perform 

the  integration  indicated  in  equations  3.13  we  proceed  by  expanding 

the  field  E   in  a  fourior  series.  We  choose  a  series  appropriate  for 
rt  r 

the  fields  over  the  reference  sphere,  in  this  case  a  scries  with  even 
symmetry  about  G'  -  0,  and  where  the  $  dependence  matches  that  cf  the 
fields  at  the  boundary,  as  suggested  by  equation  2.52.  Proceeding: 

CO 

E  ,  (R',0',<J>)  -    E    (  c  coscj)  in1   +  d  sin6  ijeosra'e?   (3.14) 
rt    *   "        rt     is      8      m      4> 
m=0  r 

The  fourier  coefficients  arc  obtained  in  the  normal  fashion,  and  arc 
covered  in  appendix  2. 

Consistent  with  the  notation  of  Ricardi  we  define  coupling  coefficients 


IC   -  /^cosmG'   P1(cosG')  dG 
ran    o         n 


(3.15) 

Id   -  /\.osr,G'  sine'   d   Pn<cos8'>   d9« 
mH    .  ____ 
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Substitution  of  3.15  into  3.13,  using  the  definition  of  fourier  co- 
efficients from  3.14  we  obtain: 


a   =   J   n       £   (  c   I    -dl) 
n     •——■*,  m  mn      m  mn 


H  (kR1)    i:  0 
n 


b   =   ~  An  R'     I      (  c   IC   -   d   Id   ) 

n     ■■  '  ■  -■■»'■■'"■•         m  mn      m  ran 

H  (kR1)   m=0 
n 


(3.16) 


where 


A 
n 


2n  +  1 
2(n(n+l))2 


Therefore  to  fully  express  the  electric  field  in  the  spherical  harmonic 
expansion  we  need  only  to  evaluate  the  fourier  coefficients  in  3.14  and 
the  coupling  coefficients  given  by  3.15.  These  coefficients,  when  sub- 
stituted into  the  scries  3. 16  yield  coefficients  of  the  spherical 

harmonic  expansion. 
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These  constants  arc  normally  evaluated  by  a  recursion  technique 

rather  than  solving  for  each  term  by  evaluation  of  the  integrals.  This 
recursion  technique  is  particularly  easy  to  realize  with  machine  computa- 
tion, and  is  inherent  in  Ricardi's  method  of  evaluation.  Thus  coefficients 
will  be  computed  rather  thaa  analytically  determined* 

1\ 


In  summary  the  relationships  that  express  the  electric  field  as  a 
spherical  expansion  of  a  specified  field  on  a  given  reference  sphere 
such  as  shown  in  figure  12  are  listed: 


,./  1    „i     ,\    -1    ^    l"*"a  n(n-'rl) 


.1 


H  (kr')   P  (cose1)   i  , 
n       n         r 


r1   n=l 


jwer' 


(3.17) 


[M2"> 
b  H  (kr' 
n  n 


A«l  1 


)  C"   +   jc  a,.  H  (kr») 


In 


n  n 


'in 


o       e 
where  C,   and  b\   are  given  by  equation  3,11  with  m-1  : 
In      In 


/n(n+l) 


P  (cose')  iQ, 


.  ,  dP  (cosG')  . 

SlU(|i   Jl  1 

dV 


3 


n(n+l) 


.  dP  (cose1)   . 

de'       e 


1 

sinji  Px(cosG')   i 


sinG' 


and  a  and  b  given  by  equations  3.16: 
n      n°      J      ^ 


1  A  R     _  ,           Td  ,   _c   x 

a  -   J   n       I  (  c   1  -  d  1   ) 

n    — — ~.~~~~.~_        m  ran  m  cm 

CHr(kr')   ra-0 


7Z 


CO 


,  -  A     II'  ,  c  .      Td      v 

b     =  n  E(cl-dl) 

H    (hr1)  n-0 

n 


2n  +  1 

A     - 


2(n(n+l))2 


c  d 

and  the  coupling  coefficients   1,1        given  by     3.15: 
r  mn'      mn 


IC     -  /cosine'    P1(cose')    dO' 
mn       o  n 


Td  ,tt        ,,    .    .,      dP    (cos9*) 

I        =   /cosm9    smO  u  d0 

mn        o  — — -- 

d0' 


The   fouricr  expansion  coef f icents   c     and  d     are  defined   in  equation   3,14: 
r  mm  ^ 


E  ,    -  E      (R1, 0'. <b)   -        E        (   c      e©sd»   i.,      +     d     sir.cb   ijcocme' 
rt         rt        '      ,y  in         Y     0'  m  <J> 

m=0 


7S 


E   is  the  vector  representation  of  the  electric  field  tangential  to  the 
rt 

reference  sphere,  as  determined  from  2,52, 

While  the  series  representation  gives  an  accurate  solution  for  the 
sum  of  all  terms,  it  is  recognised  that  this  is  impractical  to  realise. 
An  important  section  of  Ricardi's  thesis  dealt  with  justification  of  an 
appropriate  truncation  of  the  series  solution  ,  That  analysis  deals  with 
subjects  beyond  the  scope  of  the  present  work,  although  it  will  be  noted 
that  the  series  truncation  stems  from  a  desire  to  avoid  generation  of 
highly  reactive  waves  that  do  not,  in  practice,  yield  that  Much  accuracy 

to  the  end  result.  On  the  basis  of  a  rigorous  analysis  regaeding  genera- 
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tion  of  these  highly  reactive  waves  ,  it  will  suffice  here  to  state  that 

satisfactory  results  are  obtained  with  the  spherical  expansion  when  the 
number  of  terns  of  the  series  is  on  the  order  of  N  =  kR]  ,  x^hcre  R'  is 
the  particular  radius  at  which  evaluation  is  desired. 

Therefore,  equation  3.17,  evaluated  at  r'  =  R'  ,  where  R'  is  the  radius 
enclosing  the  transmitting  feed  which  will  produce  the  desired  field 


reference 
sphere 


specified  by  equation  2.52,  which  will  in  turn  produce  after  reflection 
a  desired  field  distribution  at  the  aperture,  yields 


Efeod   -   V  »£■  6''  ♦> 
sphere 
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The  equations  for  magnetic  fields  over  the  feed  sphere  can  be  obtaii 
from  Maxwell's  equations.  In  particular: 

H-   , _  -   -   *     V  x  Er       , 
feed         ~~ ——-       feed 

sphere        jtou         sphere 

In  this  analysis  vc  have,  concentrated  on  solution  of  the  electric  fields 

at  various  boundaries.  To  produce  the  desired  fields  at  the  aperture  it 

is  not  necessary  to  separately  excite  both  E..   ,   and  H..   ,   .  since 

t  feed        feed   * 

sphere      sphere 
through  Maxwell's  equations,  the  existance  of  one  field  specifies  the 
other. 
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3,2       To  illustrate  tire,  principles  outlined  in  precee.di.ng 
se.ctioi*s  and  to  show  that  tlie.se,  principles  provide  an 
economical  and  efficient  method  of  mathematical  expansion 
of  fields  from  the  geometrical  optics  reference  sphere  to 
an  interior  feed  sphere,  a  particular  case  was  chosen  for 
evaluation, 

The  parameters  of  interest  were  chosen  for  c  =  .52, 
representing  a  displacement  of  the  origin  of  the  feed  and 
reference  spheres  to  z    =  .52R  for  a  reflector  of  radius 
R  =  93A  and  subtended  angle  of  37,4  degrees.  The  field 
distribution  of  magnitude  and  phase  were  presented  in 
graphical  form  in  section  2.5  for  this  choice  of  parameters . 
Figure  3  of  section  2.1  illustrates  the  geometry  of  the 
eccentric  coordinate  system. 

It  must  be  noted  that  the  choice  of  eccentric  origin 
chosen  for  this  case  of  investigation  was  rather  arbitrary, 
as  was  the  choice  of  inner  sphere  radius.  This  choice  was 
guided,  however,  from  principles  of  geometrical  optics.  The 
rays  incident  upon  the  extreme  edge  of  the  reflector  surface 
are  reflected  into  the  focal  region,  forming  the  caustic 
surface.  The  points  at  which  the  edge  rays  are  tangent  to 
the  caustic  surface  define  the  marginal  surface  focus,  that 
plane  surface  generated  by  the  rotation  of  B-B'  of  figure 
15  about  the  axis  of  revolution  of  the  reflector. 
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Note  that  all  reflected  rays  must  pass  through  this  surface.  The 
smallest  such  surface  through  which  all  reflected  reiys  must  pass  is 
defined  by  the  circle  of  least  confusion,  generated  by  rotation,  of 
A-A1  of  figure  15  about  the  axis  of  symmetry. 

While  geometrical  optics  wrongly  predicts  infinite  fields  in  the 

vicinity  of  caustics,  the  fact  remains  that  converging  rays  indicate 

an  area  of  energy  concentration.  Hence,  an  intuitive  choice  of  probable 
f 

location  for  a  minimum-sized  feed  would  be  that  aperture  which  most 
closely  coincides  with  the  circle  of  least  confusion  associated  with 
a  particular  reflector.  It  is  a  simple  matter  to  show  from  the  geometry 
that  a  choice  of  origin  between  z  =  ,5R  and  ,6R  requires  a  feed  sphere 
on  the  order  of  5-10  X  for  a  reflector  of  93  X  with  subtended  angle 
37.4   ,  as  considered  in  this  study,  to  maet  these  specifications. 
Hence,  a  representative  distribution  of  fields  over  such  a  feed  sphere 
that  would  be  required  by  the,  tangential  boundary  conditions  to  yield 
a  uniform  plane  wave  over  the  reflector  aperture  is  given  by  figures 
16.   Amplitude  and  phase  of  three  components  of  electric  field  over 
the  surface  of  a  feed  sphere  of  radius  8.4  X  with  origin  at  z  =  ,52R 
for  a  reflector  of  radius  R  ~  93  X  are  presented  in  figures  16  A-F. 
This  solution  is  a  graphical  result,  obtained  from  evaluation  of  the 
series  solution  given  by  equation  3.16  and  3.17  of  section  3.1. 
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3,3   For  the  choice  of  parameters  indicated  in  section  3,2 
the  magnitude  and  phase  of  the  field  distribution  over  the 
reference  sphere,  as  given  by  equation  2.5  2,  were  presented 
in  graphical  form  as  figures  13  and  14  of  section  2.5.  The 
nearly  identical  distribution  of  tangential  components  and 
the  near ^vanishing  radial  component  of  the  field  are  due 
to  the  relationship  between  angular  variables  8  and  0', 
where  to  first  order  ; 
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Only  the  tangential  field  components  are  used  in  the 
boundary  condition  to  generate  the  fourier  coefficients, 
as  noted  in  appendix  II.  The  tangential  fields  were  expanded 
in  the  following  manner  from  equation  2.52: 


E  =     (       E    cos<f>i     ,    +    E    sin<J>i       )ej4\ 

tan  0099 


E^.  =    E„  cosrf>  (cosiL    +    jsinij;)     in    +    E.sin^cosijj    +    jsini|;)i 

tan  0  0  9  9 
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The.  real,  a.nd  imaginary  pa,rts,  le.ss  th^e  cos4>?  sin<f>  dependence 
were  then  plotted  graphically  and  are  presented  as  figures 
17  and  18  ,  Due  to  the  similarity  of  magnitudes  resulting 
from  the  nearly-identical  nature  of  tangential  components 
indicated  by  figures  13  and  14  of  section  2.5,  only  two 
rather  than  four  plots  are  presented  here.  The  calculations, 
however,  retain  the  existence  of  four  separate  components, 
two  real  and  two  imaginary,  of  tangential  fields.  These 
components  were  then  expanded  into  a  fourier  series  of  the 
form : 


M 

E  ^  =       ^n(    c       cos|    i0,    +    d       s  i  n  <j>    i,)cosmO 

tan         m=0 v      m  6  '  m  <J> 


where  c  ,  d   are  now  complex  fourier  coefficients  determined 
mm 

as  noted  in  appendix  II,  and  M  is  the  limiting  value  of 
the  highest  harmonic  desired,  given,  as  noted  previously, 
by 


M 


kr 


with  r  the  radius  of  observation 


To  check  on  the  accuracy  of  the  fourier  expansion  technique,  the  series 

was  summed  using  coefficients  determined  by  the  geometrical  optics 

field  over  the  reference  sphere.  The  results  are  indicated  in  figures 

17  and  18  as  small  departures  from  the  heavy  dark  line,  which  represents 

the  values  used  to  generate  the  fourier  coefficients. 

These  fourier  coefficients  were  then  used  vzith  the  series  solution, 

given  by  equations  3.16  and  3.17,  to  evaluate  the  fields  at  the  reference 

sphere.  The  trends  and  accuracy  of  this  series  solution,  when  compared 

to  the  geometrical  optics  field  initially  calculated  and  subsequently 

used  to  generate  the  necessary  coefficients,  appears  to  be  reasonable 

when  viewed  in  figures  19.  It  is  to  be  recalled  that  the  series  was 

truncated  at  a  value  of  terms  much  Decs  than  infinite  in  order  to  avoid 

generation  of  highly  reactive  fields.  The  presentation  of  figures  19 

involved  390  terms.  As  the  series'  accuracy  improves  with  the  number  of 

terms,  better  agreement  between  the  geometrical  optics  and  resulting 

synthesized  fields  would  theoretically  be  forthcoming  with  continued 

summation.  It  is  to  be  noticed  that  the  agreement  of  phase  appears 

quite  good  for  the  tangential  fields  and  satisfactory  for  the  radial 

component.  Similar  agreement  occurs  in  amplitude.  The  basic  boundary 

condition  that  was  involved  required  matching  of  the  tangential  fields 

so  as  to  force  the  field  to  zero  at  tlie  reflector  surface.  Hence,  the 

radial  component  is  also  determined  from  this  boundary  condition,  and 

should  contain  some  variation  for  less  than  an  infinite  number  of  terms. 

The  phase  of  the  radial  component  is  similarly  determined  from  the 

real  and  imaginary  parts  of  the  component,  and  for  the  radial  component 

89 


of  near-zero  magnitude,  the  error  involved  in  the  difference  of  small 
numbers  begins  to  override  the  solution,  Hence,  for  values  of  radial 
component  near  zero,  rapid  phase  variations  were  encountered.  These 
were  not  plotted  as  not  being  significant  they  added  no  contribution 
to  the  work. 
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CHAPTER  IV 
CONCLUSIONS 


4.1  As  noted  by  the  graphical  representation  of  solution  offered  in 

figures  19,  the  technique  explored  by  this  thesis  has  yielded  results 
that  indicate  its  accuracy  and  applicability.  A  geometrical  optics  expan- 
sion from  the  surface  of  a  spherical  reflector  to  a  hypothetical  ref- 
erence sphere  can  be  utilizedto  generate  necessary  coefficients  required 
by  the  series  solution  of  spherical  harmonics  that  constitute  the 
expression  for  the  electromagnetic  fields  in  the  region  of  concern. 

With  regards  to  the  efficiency  of  the  technique,  an  investigation  of 
various  eccentric  origins  from  z  -  ,5R  to  .6R  indicated  that  the  radii 
of  the  ensuing  reference  spheres  were  on  the  order  of  R'  ~  ,6R  for 
the  reflector  chosen  for  investigation  by  this  thesis.  Thus,  evaluation 
of  the  series  solution  at  the  reference  sphere  for  purposes  of  compar- 
ison to  the  initially  determined  geometrical  optics  field  involved  a 
reduction  in  the  number  of  terms  that  would  have  been  required  without 
the  eccentric  coordinate  system  by  some  40  %.  The  reduction  in  terms  is 
due  to  the  direct  relationship  between  the  radius  of  observation  and  the 
number  of  terms  required.  Due  to  the  flexible  nature  of  the  technique, 
choice  of  eccentric  origin  to  include  regions  of  energy  concentration 
near  the  feed  sphere  result  in  the  feed  radius  on  the  order  of  8  X 
for  the  reflector  of  93  X   considered  by  this  study. 
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Previous  investigation    constrained  by  a  fixed  origin  and  con- 
centric spheres,  as  given  by  figure  2,  found  it  necessary  to  keep  the 
radius  of  the  feed  sphere  on  the  order  of  half  the  radius  of  the  reflec- 
tor* Thus  we  consider  that  this  technique  has  provided  a  substantial 
reduction  in  the  number  of  terms  involved  from 


M  =  kr  »  (  2-.t/X  )(  93X  /?.) 


to 


M  -  kr  =  (  2t./X  )  8  X 


indicating  a  savings  of  a  factor  of  approximately  an  order  of  magnitude. 
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4.2   A  next  1  o g i  c a  1  ste.p  in  p u r  s u a n c e  of  t hi  s  tec hn ique 
would  be  to  concentrate  on  the  series  solution  at  the 
radii  of  various  feed  spheres.  This  study  suggests  that 
if  the  fields  given  by  figure  16  can  be  duplicated  over 
the  surface  of  the  feed  sphere  by  some  source  within  the. 
sphere  the  desired  illumination  of  uniform  plane  wave 
over  the  aperture  of  the  reflector  will  be  realized.  The 
calculation  of  the  fields  at  the  aperture  by  an  alternative 
measurement  technique,  or  by  experimental  methods,  if  the 
fields  over  the.  feed  sphere  can  be  realized  by  some  practi- 
cal feed  system,  would  be  a  most  worthwhile  endeavor. 

Further  study  for  a  series  of  parameter  variations  for 
one  reflector  could,  as  well,  yield  a  meaningful  extension 
of  this  technique.  Based  upon  the  field  distribution  over 
the  feed  sphere,  design  and  experimentation  with  a  feed 
system  could  provide  the  basis  for  additional  study  of  this 
topic  at  various  academic  levels  of  concern. 


33 


APPENDIX  I 
SIMPLIFICATION  AND  MANIPULATIVE  DETAILS  OF  CERTAIN  EXPRESSIONS 


D a  r  iv a  t  :i  o  n  o  f  Eo  ua t  ion  2  c  1 2 


F  r  c  m  c  q u a t  i  o n  2.11: 


N/p  =  -  S  x  7  x  S  A.l 


Taking  the  dot  product  with  N  to  both  sides  of  A.l 


1/p  =  ~N*S  x  v  x  S 


applying  the  vector  identity 


A«(B  x  C)  =  C(A  x  B) 


then 


1/p  -  -  (N  x  S).(7  x  S)  A. 2 


By  equation  2. 10 


S  -  7L/n 
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Substituting  and  operating: 


V  x  S  -  V  x        -  V  1        x  VL  -!•   -1   V  x  VL 


But   V  x  VL  -  0  so 


V  x  S  -  V   1   x  VL  •        A. 3 


Substituting  A. 3  into  A, 2: 


1/p  =  -(N  x  S).(  V  l        x  VL)  A. A 


We  substitute 


v  _L-   =   ~    X   V  (In  f() 


1/p  «  -(N  x  S).(  -  J_   V(ln  f,  )  x  VL)       A. 5 

n 


But  by  the  vector  identity 


A  x  B  •  C  x  D  -   (A«C)(B»D)  -  (A«D) (B«C) 


»0| 


Ac6  1/p  -    1  (  (N-Vln  n)(S-7L)  -  (N»7L)(S-71n  n)  ) 


Fram  2.10,  where  S*VL  -   |vl[2  «  n 

n 


1/p  «      (  N-Vln  n  -  (N«VL)'(S»Vln  n)      A. 7 


The  second  tern  in  A. 7  is  identically  zero  since  by  A«l 


N  =  -  p(S  x  V  x  S) 


and  fron  2.10 


7L  *=   Sn 


then  substituting  the  above  into  the  second  term; 


(N«7L)  =  -S«(S  x  V  x  S) 

pn 


Again  using  the  vector  identity 


A»(3  x  C)  «  C«(A  x  B) 
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we  GGC 

V  x  S  «(S  x  S)  «  0 

Therefore,  from  Ac7 


1/p  =  N»Vln  n  A. 8  »  2.12 


Derivation  of  equation  2,37  ' 


From  equation  2,34 


x  =  (  z  -  z   )  tan20 


x2=  (  z2  -  2zz  +  z2  )  tan2 29  A. 9 

o    o 


From  equation  2,35 


R'2  -  x2  +  z2  -  2zcR  +  c2R2  A. 10 


Substitution  of  A. 9  into  A, 10 


R'2  =  (  z2  -  2zz  +  z2   )tan220  +  z2  -  2zcR  +  c~i 


2^2 
o    o 


A,  11  =(tan226  +  l)z2  •-  2z(z  tan228  -I-  cR)+  z2tan220+c2R2 
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Multiplying  A. 11  by  cos229 


Rf2cos220  -   (sin229  +  ccs229)z2  -  2z(z  cin220  +  cRcgg229) 

© 


+  z2  sin22B  •■!-  c2R2ces226  A.  12 

o 


Solving  A. 12  for  z  i 


0  -  z2   -   2z(z  sin220  +  cRcos229)  +  z  2sin229  -  cos226(R*2  -  c2R2)   A. 13 
o  o 


From  2.36; 


R  z 

sin29         sine 


sine  _,  .  , , 

z      =  R  A.  14 

Din  26 


Substituting  A. 14  into  A. 13  : 

z2  -  2z(sin9sin20R  +  cRcos226)  +  sin29R2  -  cos229(  R'2  -  c2R2)  =  0 

A, 15 


wi 


Normalizing  A. 15  to  the  radius  of  reflector  R: 


Z  n.        & 

z     ~  t      R'  - 

n   — -   f   n    — 


.  t 


R 


R 


z2  -  2z  (sin8sin26  +  c2coc22G)  +  sin26  ~cos229  (R'2  -  c2)  -  09A. 16=2.37 
n     n  n 


Derivation  of  equations  2 . 43 ,  2.44 


Normalizing  equation  2.41,  2.42: 


R,      =   /(   z       -  z     J2"*   (  x       -  x     )• 

In  mn  en  ran  en 


A. 17 


v?here 


z        =   co;; 
ran 


x       =  sine 

inn 


en 


1_(   3cos0   -   cos30) 
4  ' 

x        =  (   3sinG   -   c.in36) 

en       —  7  ■— 
4 


z        =      sinO 
on  .    -  ' 

sin20 


Subs  t i tut i on     and  cxp  an d i a g  A . 1 7 : 


/OS 


R2        «      (     cesO  -!-  cr   ?><■:■   )?   ■:-      (     ^1:0     -:   ;  '    ' '■';    )2  A. 18 

1^  *   '/,  A  "    /.  A 


-   1    (  2c©s8c©s3e  +  2sin9sir.30  +  2  ) 
16 


2^_  (  ?  +  2cos20  )  A.  19 

16 


but   cos2G  "  1  -  2s in' 


Ri2  "  ^Jl^S   1  "  siu20) 
Jn      16 


R2   =  cos2e 

in     ——; — — 


then 


Ri    =  -^2£L-  A.  20=2.41 


Similarly,  by  expanding  2.42: 


R2    =   (  cos0  _   sinG   )2  +  sin2e  A%21 

211  7in?9 


~      cos?0  -!-  sin2  6  +   sd       -  _  "/ 

sin2 2 9        sin26 


lo& 


R 


2  n 


1   -   1   +     _sj 


y 


R 


2n 


A. 22-   2.44 


2   con9 
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APPENDIX  IT- 
COMPUTATIONAL  TECHNIQUES  RELEVANT  TO  PROBLEM  SOLUTION 


As  noted  in  section  2.5s,  expressions  for  fields  were  programmed  in 
FORTRAN  to  yield  a  numerical  solution  to  equation  2,52,  A  copy  cf  the 
program  utilized  follows  this  discussion.  Briefly,  the  logical  progress ion 

is  outlined: 

1.  For  a  given  reflector  of  radius  R  and   subtended  angle 
On,  equation  2,31  is  used  to  calculate  the  reference 
sphere  radius  R1  for  a  particular  choice  of  parameter 
c, 

2.  Equation  2.40  is  programmed  to  yield  the  angle  0' 
ueasured  from  the  new  origin  at  2=cR  to  points  on  the 
reference  sphere  corresponding  to  rays  incident  upon 
the  reflector  surface  at  a  given  angle  0,  A  one-to-one 
correspondence  is  established  between  the  angles, 

3.  Equations  2.61  and  2.62  are  programmed  along  with 
2,33  to  yield  values  for  the  field  amplitude  at  var- 
ious angles  G',  measured  from  the  new  origin,  en  the 
surface  of  the  reference  sphere.  The  ai  plitude  expression 
is  then  given  by  the  inverse  distance  relationship  of 
the  form  of  equation  2,17. 


A.  Th«  polarization  features  ccvcrc,'.  by  equations  2C/;6"2.51 
are  prog.,     I  to  indicate  the  effect  of  redefining  the 
field  cc;  j    ats  in  tl      e  »vdi   tate  syst<  ■ 

5«  The  phase  expr  i  ion  is  calculated  for  the  reflector 
under  consideration  by  relation  2.30  and  2.33,  and  V 
resulting  phase  is  chosen  relative  to  the  pi     at 

0  -  0f  ^  0  degrees. 

6.  The  information  generated  by  the  previous  steps  is 
collected  and  output,  if  desired,  is  Hade  available, 

representing  a  numerical  solution  to  equation  2,52. 

The  graphs  included  in  section  2.5  were  prepared  from  results  generated 
as  described  above. 

It  is  noted,  however,  that  the  actual  field  on  the  reference  sphere 
is  only  of  passing  interest  at  this  point,  and  that  expansion  of  this 
field  into  a  cosine  fourier  series,  as  given  by  equation  3.14,  is  cf 
paramount  interest: 


E  reference  =   E   (  c  cos<j>  i. ,  +  d  sin*  i^ccsp.9     A. 23 
sphere         ra=0  Y 

tangential 
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The  subroutine  FORC  was  e!  .  h   pad  to  generate  coefficients  for  the 
fourier  series  front  the  generated  field  c  .      ,  A  numerical  technique 

for  evaluating  the  coefficients  for  rich  a  series  appropriate  for 

A  A 
a  periodic  set  cf  discrete  values  of  a  fu       was  available  r.\u\   was 

used  as  the  basis  for  this  subroutine.  As  a  check,  the  fourier  series 

was  si;; .  .e'.1,  i:  ■■■>.(:.      '  g  accuracy  to  within  a  few  percent  when  viewed 

overall  and  compared  to  the  points  used  to  generate  the  coefficients. 

The  technique  follows  fro^a  consideration  of  expansion  of  a  periodic 

function  f(0)  in  a  fourier  series; 


f  (6)  -     a  /2  +  I   a  cesm6  +  b  sir..  18  A«24 

o      .  ta  t.\ 


where  from  ertheg   lity  it  nay  easily  be  shown  that  the  resultant 
coefficients  arc  given  by: 


2tt 
I 
0 


a  -        /  f(G)  cesrae  d.6  n'-O,  19  ....     A, 25 

la     — — = 


As  the  case  of  interest  involves  symmetry  about  0-0,  it  is  obvious 

that  a  cosine  series  will  suffice  to  fully  express  the  function  and 

b  will  be  identically  zero  for  all  v.. 

HO 


The  integral  expression  A, 25  can  be  related  to  a  su  tatit  a  of  discrete 
values  by  app:  .  Lmatiens  of  basic  calculus,  where  the  interval  fr< 
0  <  0  <  271  is  divided  into  2N  intervals,  providing  2N  -!•  1  data  points, 
Thus  we  can  formulate  a  relationship  between  the  number  of  intervals 
and  the  spacing  involved  where 


2N  =   360°/  DEL 


DEL  being  the  spacing  between  successive  points .Thus  we  approximate, 

for  a  sufficiently  large  number  of  points: 


dO  *  A6  «   27T 

2N  +  1 


where  n~-  0,1,2,  ....2:1 


2N  +  1 


Substitution  of  the  above  approximations  into  A, 25  results  in 


£n  g  £    *<_„>  eosm   """ 

2N  +  1 


IK 


Thus  we  ceo  the  function  must  be  na.de.   available  in  steps  of   angle, 
determined  by  the  index  n,  or  at  the  following  discrete  points: 


S   -     2irn        n  -  0,1,2,.  ,.2N 

2N  +  1 


Therefore 8  the  following  steps  are  outlined  to  provide  an  indication 
of  the  progression  followed  through  the  subroutine: 

1«  Since  the  scries  requires  even  spacing  between  the  data 

points,  it  was  necessary  to  interpolate  between  tbe  points- 
generated  on  the  reference  sphere  by  successive  steps  of 
the  angle  along  the  reflector  surface.  As  equation  2.40 
allowed  no  sinple  inverse  relationship  between  the  two 
angles  0  and  0',  snail  steps  of  the  angle  on  the  reflector 
surface  were  taken  (  spacing  -    .2  degrees  )  to  insure  the 
. accuracy  of  linear  interpolation  between  the  discrete  values 
of  angle  0'  generated. 

2.  The  resulting  field  expressions  were  arranged  in  component 
fom  and  subsequently  expanded  fro:i  the  exponential  represent- 
ation to  the  trigonometric  fom  : 


HZ 


where  Ae   -  A(  cos^  +  jsinip  )  -  AB  +  jAI 


3,'  The  values  of  field  c-.  ;-.  ei  ts,  twe  real  and  two 
imaginary,  were  then  used  ,te  g     :   a  c  plex 
fcurier  coefficient  of  the  fore 

C  -  CR  +  jCI 
by  the  technique  previously  described. 

4*  The  answers  of  interest  were  these  coefficients  \</hich, 
when  properly  inserted  into  an  existing  program  for 
the  series  solution  to  field  expressions  3.16  and  3.17. 
The  fourier  coefficients  c  and  d  ere  noted  in  equation 
3.16. 
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